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The Haldane pseudopotential construction has been an extremely powerful concept in quantum 
Hall physics — it not only gives a minimal description of the space of Hamiltonians but also suggests 
special model Hamiltonians (those where certain pseudopotential are set to zero) that may have 
exactly solvable ground states with interesting properties. The purpose of this paper is to generalize 
the pseudopotential construction to situations where interactions are A-body and where the particles 
may have internal degrees of freedom such as spin or valley index. Assuming a rotationally invariant 
Hamiltonian, the essence of the problem is to obtain a full basis of wavefunctions for N particles 
with fixed relative angular momentum L. This basis decomposes into representations of SU(n) with 
n the number of internal degrees of freedom. We give special attention to the case where the internal 
degree of freedom has n — 2 states, which encompasses the important cases of spin-1/2 particles 
and quantum Hall bilayers. We also discuss in some detail the cases of spin-1 particles (n = 3) and 
graphene (n = 4, including two spin and two valley degrees of freedom). 

PACS numbers: 73.43.-f, 03.65.Fd. 



I. INTRODUCTION 

In the absence of interactions or disorder, the impact of 
a magnetic field on the band structure of two dimensional 
electrons is profound: the spectrum of single particle 
cigenstates breaks into degenerate bands called Landau 
levels (LLs). Fractional quantum Hall (FQH) physics^ 
occurs when interactions between electrons break the de- 
generacy of partially filled Landau bands leading to an 
incompressible fluid ground state. 

In most discussions of fractional quantum Hall effect 
(FQHE), one further assumes that the interaction be- 
tween particles is a pair- wise potential V(|ri — r^l) and 
that the electrons have no additional degrees of freedom 
such as a spin or valley index. These two further assump- 
tions are ones that we would like relax in the current 
work. In this paper, we study so-called multi-component 
systems (particles with an additional degree of freedom, 
such as spin) interacting via very general iV-body poten- 
tials V(ri, si; r?, S2', ■ ■ ■ rjv, sjv) where here Sj represents 
the additional quantum number of the j particle. Mo- 
tivation for expanding our attention to this more compli- 
cated case is discussed below in section ITBl 

In analysing FQH physics a very important tool is 
the Haldano^ pseudopotential. In the simple case where 
the interaction is two-body and there is a single species 
of spin-polarized electrons (i.e., without the generaliza- 
tions discussed in this paper) , the Haldane pseudopoten- 
tial, Vl, is defined as the energy cost for two electrons 
to have a given relative angular momentum L. These 
pseudopotentials give a complete (and minimal) descrip- 
tion of any rotationally and translationally invariant two- 
body interaction within a single Landau level. The lan- 
guage of pseudopotentials not only provides a convenient 
parametrization of the problem, but also makes it easy to 
write down certain model Hamiltonians which are solv- 
able — thus providing a key piece of our understanding 
of the FQHE. The purpose of the current paper is to 



generalize the idea of the pseudopotential to the multi- 
component case where particles have additional degrees 
of freedom such as spin or valley index, and where inter- 
actions between particles are of a general A^-body form. 



A. Pseudopotentials and Quantum Hall States 

Let us elaborate for a moment on the connection be- 
tween pseudopotentials and exactly solvable model inter- 
actions in the case where interactions are of the conven- 
tional two-body form. Let us choose a model interaction 
potential such that the value of V\ is positive, but Vl — 
for L > 1. In this case, we are imposing an energy cost 
for any two particles to have relative angular momentum 
of L = 1 : the exact highest density zero energy state (the 
highest density state where no two electrons have relative 
angular momentum of 1) is precisely the v = 1/3 Laugh- 
lin wavefunction^— (as always, we define filling fraction 
v to be the fraction of single particle states within a LL 
which are occupied — or equivalently v is the ratio of 
number of electrons to number of flux quanta). Because 
of this exact solvability, much can be established in de- 
tail about this type of model Hamiltonian. While the 
model Hamiltonian may not be too similar to any partic- 
ular physical Hamiltonian (such as the Coulomb Hamil- 
tonian), nonetheless, the ground state may be very close 
to the physical one. More importantly, the ground state 
of the model gives an easily studied representative of an 
entire phase of matter. 

In recent years, due to a large extent to interest in 
more exotic non-Abelian quantum Hall states^, interest 
has turned from two body interactions (such as Coulomb) 
to A^-body interactions with N > 2. While in principle 
such Af-body interactions occur as a result of integrat- 
ing out higher Landau levels^, and could in principle be 
engineered to exist in certain cold-atom systems^, the 
main interest in multi-body interactions is again due to 
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the fact that certain interesting wavcfunctions are the 
exact highest density zero energy state of certain N- 
body interactions. For example, the celebrated Moore- 
Read wavefunction 8 - is the exact ground state of a model 
three-body interaction 9 -, and more generally the Read- 
Rezayi wavefunctioniS is the ground state of a (k + 1)- 
body interaction. Other interesting examples include the 
Gaffnianii and Haffnian^ wavcfunctions which are also 
ground states of special three-body interactions. 

As with the simple case of two-body interactions, 
when studying FQHE Hamiltonians with TV-body inter- 
actions, it is quite useful to work with a generalization 
of Haldane's pseudopotential formalism first proposed in 
Ref. [TH While this formalism is similar to that dis- 
cussed by Haldane in that it decomposes the interaction 
into angular momentum components, it is a bit more 
complicated since specifying the total relative angular 
momentum of N particles does not completely specify 
the "relative" wavefunction of the iV-particlcs as it does 
in the case of two particles. This complication is ex- 
plained in detail in section |TT] below where we review 
the multi-particle pseudopotential formalism. Account- 
ing for this complication, the multi-particle pseudopoten- 
tials can simply describe all of the special model Hamil- 
tonians pointed out above. In fact, in each of the above 
discussed cases, the Hamiltonian can always be described 
as simply forbidding any N particles from having relative 
angular less than some L. 

All of the wavcfunctions discussed above turn out to 
be special cases of a very broad set of so-called Jack 
polynomial wavefunctions^. While there is a general be- 
lief that local ./V-body Hamiltonians may exist for all of 
the Jacks, only limited further cases have actually been 
explore d 15 ! 16 . Such Hamiltonians, should they be con- 
structed, can be phrased in the general language of pseu- 
dopotcntials. In these more general cases, the necessary 
Hamiltonians cannot be of a simple form that simply pro- 
hibits clusters of particles from having certain angular 
momenta^, but rather will require a particular balance 
of certain pseudopotential coefficients^. 



B. Multi-Component Quantum Hall 

All of the FQH states discussed so far describe spin- 
polarized or "spinless" system — i.e., the only degree of 
freedom for each electron is its orbital position. However, 
there are many cases where we may want to consider a 
more complex model where electrons are endowed with 
additional degrees of freedom, such as spin, and the FQH 
ground state may not be fully polarized. In these situa- 
tions we say that the ground state is "multi-component" , 
meaning that it involves an important contribution from 
electrons having different values of these additional in- 
ternal degrees of freedom. 

The simplest example of a multi-component quantum 
Hall system is one with spin. Naively, one might expect 
that the high magnetic field characteristic of the frac- 



tional quantum Hall regime might remove any spin degree 
of freedom entirely. However, in many quantum Hall sys- 
tems this is not the case, and we must consider the spin 
degree of freedom as well. For example, in conventional 
GaAs systems, owing to the small g-factor, even in fairly 
high fields, electron spin may not be polarize d 18 ' 19 . Fur- 
thermore, there are experimental methods^ to engineer 
an even smaller g-factor, making both spin states even 
more relevant. 

Valley degrees of freedom are another way which quan- 
tum Hall systems may be multi-component. For exam- 
ple, in both AlAs quantum wells^i, and graphene^r— 
the semiconductor band structure is such that each elec- 
tron has a spin and valley degree of freedom. Analogous 
to the electron spin, in these cases, the valley index may 
take one of two possible values, so we think of the valley 
as a "pseudospin" or "isospin". Silicon MOSFETs may 
be even more complicated — depending on the crystal 
orientation with respect the the 2d electron layer, the 
electrons will have a valley index which takes one of 2, 4, 
or even 6 values 2 ^—. Yet another important case where 
electrons have an additional degree of freedom is in quan- 
tum Hall bilayers, where the layer index plays the role of a 
pseudospi n 28 ' 29 . All of these multi-component quantum 
Hall systems have been the subjects of intense theoretical 
and experimental study. 

A related topic of recent theoretical interest is the 
study of the quantum Hall effect of cold bosons. While 
such a quantum Hall effect has not yet been observed, 
it seems potentially feasible in the not-so-far future 3 ^. 
Bosons, as compared to electrons, must have integer 
spin, but one could easily imagine a quantum Hall effect 
of spin-1 bosons which would have three internal states 
rather than the two of an electron. Other possibilities for 
multi-component Bose systems exploit multiple hyper- 
fine states of an atom, or multiple sub-bands that occur 
for bosons in a magnetic field on an optical lattic o 31 ' 32 . 
Another application would be to systems where multiple 
Landau levels can be occupied where the internal degree 
of freedom would be the Landau level index. 

The many possibilities of experimentally realizing 
multi-component quantum Hall states have driven a large 
number of theoretical studies^ 8 -. These studies began 
with very early work by Halperin 3 ^ which generalized 
Laughlin's wavefunction^ to the multi-component case. 
Quite naturally, the concept of Haldane's pseudopoten- 
tials were quickly generalized to the multi-component 
case as well^ir— . Analogous to the Laughlin case, the 
Halpcrin wavcfunctions can also be described as the exact 
highest density zero energy states of special (two-body) 
interactions, and these interactions can in turn be de- 
scribed by particular multi-component pseudopotential 
coefficients. 

The added richness of multi-component systems has 
made them a prime place to search for new and excit- 
ing physics. In the search for novel non-abelian quantum 
Hall systems 5 -, several multi-component candidates have 
been proposed 3 ^—, including the so-called non-abelian 
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spin singlet (NASS) states 3 -^ and the spin-charge sepa- 
rated states^ for the two component case, as well as 
generalizations of these constructions to higher numbers 
of component a 42 ' 43 . Recently the idea of Jack polynomi- 
als was generalized to certain multicomponent systems 
as well 4 ^ 

Analogous to the situation with single component 
wavefunctions, many of these novel multi-component 
wavefunctions are exact ground states of special ./V-body 
interactions. Here, however, interactions may be more 
complicated — depending on the "spin-state" as well as 
the position of the particles. As such, it seems natural 
to try to generalize the pseudopotential formalism to the 
case of multi-component N body interactions. This is 
the aim of the current paper. 

As in the single component case, the motivation 
for developing the pseudopotential formalism for multi- 
component many-body-interactions is several-fold. On 
the one hand, pseudopotentials provide a complete 
paramctrization of the problem for such systems. The 
usefulness of this is evidenced by recent works which have 
introduced pseudopotentials for multi-component many- 
body interactions for important special case a 6 ' 46 . More 
importantly, the pseudopotential structure hints at what 
sorts of simple Hamiltonians may be written such that 
interesting quantum Hall states might be found as the 
highest density zero energy state of a particular set of 
pseudopotential coefficients. While in the current paper 
we do not yet undertake to identify new wavefunctions in 
this way, the research program is nonetheless clear and 
will be a topic of future research (see also Refs. HH and 

E). 



C. Structure of this Paper 

In section|TT]we review the idea of (multi-particle) pseu- 
dopotentials and we define our problem in more detail. 
Following this necessary background, our main results 
are summarized in section IIIII While our methods are 
generally applicable to any number of internal degrees of 
freedom per particle, we emphasize in particular the most 
experimentally relevant multi-component cases which are 
the cases where particles (fermions or bosons) have 2, 3, 
or 4 possible internal states for interactions between small 
numbers (2-5) of particles. We present tables for these 
simple cases indicating the number p of linearly indepen- 
dent wavefunctions that exist for N particles having total 
relative angular momentum L. This defines a p x p Her- 
mitian matrix of pseudopotential parameters that can be 
defined for that value of angular momentum. In this same 
section wc further give the explicit form of these linearly 
independent wavefunctions. For many readers who are 
interested in applications of our work, these tallies in the 
results section should provide most of the relevant infor- 
mation. Within this section, subsection IIII Al addresses 
the case of two-component wavefunctions (applicable to 
spin-1/2 systems, bilayers, etc). This section is the most 



straightforward of our results and is also probably of the 
most wide interest. Subsection IIII Bl gives a brief primer 
on symmetry types and Young tableaux which is neces- 
sary for the description of the extension of this work to 
higher number of components, which we give in subsec- 
tion IIII CI Following the results section we give a brief 
discussion in section IIVI of applications of our results as 
well as some simple examples and wc direct our readers to 
this section for a more general conclusion and discussion 
of results. 

The main mathematical formalism that derives the re- 
sults presented here is relegated to the appendices. For 
the interested reader, these appendices have been made 
fairly extensive and pedagogical. However, for most read- 
ers who are interested in the application to quantum Hall 
physics, the group-theoretical details of the derivation 
will not be necessary. Although this presentation may 
be a departure from typical structure of most publica- 
tions, wc hope that it provides the clearest approach. 

II. REVIEW OF MULTI-PARTICLE 
PSEUDOPOTENTIALS 

To pedagogically introduce the concept of the pseu- 
dopotential, we start by examining the single-component 
situation where the particles have no internal degree of 
freedom. Wc will further simplify to two-body interac- 
tion then generalize to multi-body interaction. Much of 
this exposition follows that of Ref. [H] and we refer the 
reader to that reference for more details. Finally, we will 
introduce multi-component wavefunctions in section IIIII 
below. 

Before studying interactions between particles, we 
must first discuss wavefunctions for nonintcracting parti- 
cles in a magnetic field. In the lowest Landau level (LLL) 
with a planar geometry, the nonintcracting Hamiltonian 
is 



where in the symmetric gauge A = — \v x B, with the 
magnetic flux density given by B = Bz. The solutions 
of the Schrodingcr equation in the LLL are given by 



V27r2 m m! 

where z = (x + iy) is a complex number representing 
the position of the particle in the plane, and we have 
set the magnetic length {he/ eB) 1 / 2 to unity. Here, the 
angular momentum around the origin is km and m > 
is an integer. On other geometries such as the sphere, 
solutions still will take the form z m only changing the 
normalization and the Gaussian measur o 10 ' 13 . 

More generally, a solution to the single particle 
Schrodinger equation is any linear combination of the 
above basis states, thus taking the form 

■4> = f(z)e- L ^\ 
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where f(z) is any analytic function of z. Note that due 
to an isomorphism between Landau levels it will be suf- 
ficient to study our problem in the lowest Landau level 
(LLL) only. Interacting systems within higher Landau 
levels may be treated with appropriate transformations 
of the lowest Landau level (see Ref. [ID for a detailed dis- 
cussion). We also comment that if the particle has an 
additional degree of freedom, such as spin, which we dis- 
cuss below, we will need to specify the state of this spin 
as well. For example, we might write 

4> = f(z)e-^ 2 |t> 

For the many-body version of the problem we may 
construct linear combinations of products of the single 
particle wavefunctions in the co-ordinates Zi, and im- 
pose an appropriate symmetry (a wavefunction describ- 
ing fermions must be overall fully antisymmetric, and a 
wavefunction describing bosons must be fully symmet- 
ric) . In the case of a single component wavefunction (no 
internal degree of freedom) one can simply write 

1 \— ^ N I 1 2 

V> = 0(zi,...,ZAr)e"*^=i |Zl1 

where (f> is an analytic function of all of its arguments 
and is overall symmetric for bosons and antisymmetric 
for fermions. In the multi-component case, the symmetry 
condition is much more complicated as we will see below. 
For now, we will continue to focus on the simpler single- 
component case. 

A. Single-Component Two-Body Interactions 

We begin our discussion of interactions with the case of 
simple two-body interaction. Very generally a wavefunc- 
tion for two particles can be decomposed into relative and 
center of mass components. For two-body wavefunctions 
within the LLL, we can write a complete basis 

|* L>; ) = |0 ® \L) 

where I is the center of mass angular momentum and L 
is the relative angular momentum (i.e the relative an- 
gular momentum about the common center of mass of 
the N particles). Explicitly we mean the center of mass 
wavefunction is given by 

|Z> OC + ^2)' e -|l^+^l 2 

and the relative wavefunction is 

\L) o,( Zl -z 2 ) L e-h^^\ (1) 

It is important that these wavefunctions form a complete 
set. Physically we deal with particles which are either 
fermions or bosons and so in fact the space of two particle 
eigenstates exists only for odd L for fermions (or even L 
for bosons) in order to obey the correct symmetry. 



After projecting to a single Landau level (which can 
be justified by a large gap between Landau levels) the 
Hamiltonian is simply the interaction between the parti- 
cles. We can thus write 

H = J2v(\r i -r k \). (2) 

j<k 

where V is some inter-particle interaction such as the 
Coulomb interaction. To decompose this interaction into 
pseudopotentials, we define 

V L = (L\V(\r 1 -r 2 \)\L). (3) 

Since the interaction is translationally invariant, it is in- 
dependent of the center of mass degree of freedom of the 
two particles. 

Wc can then rewrite the Hamiltonian given in equation 
© as 

H = J2J2\L)V L {L\, (4) 

L i<j 

Here it is implied that in each term of the sum, the two- 
particle ket \L) involves particles i and j. Again, since 
the interaction is translationally invariant, interaction 
between particles i and j never changes their common 
center of mass, so we may work with only the relative 
wavefunctions. Further, the interaction Eq. [3] and hence 
the Hamiltonian Eq. H]is diagonal in L on account of the 
fact that the interaction potential is rotationally invari- 
ant and therefore conserves relative angular momentum. 
The intuition behind this rewriting of the Hamiltonian is 
that any two particles with relative angular momentum 
L incur an energy cost Vl ■ 

B. Single-Component Multi-Particle Interactions 

In order to extend this discussion to multi-particle 
pseudopotentials, we consider a general many body inter- 
action potential V (ri, rjy), and we restrict ourselves 
to a rotationally and translationally invariant system. 
Analogous to Eq. [2] the LLL Hamiltonian is written in 
terms of the A-particlc potential as: 

H= ^(ri,...,rjv). 

il<»2<-"<»N 

By analogy with Haldanc's original pscudopotcntial con- 
struction, wc decompose the wavefunction of the N par- 
ticles into a center of mass and a relative motion. Further 
we would like to write a complete basis for the possible 
relative wavefunctions of the N particles which we will 
use for our pseudopotential construction. 

Determining this complete basis turns out to be the 
tricky part of the .ZV-body pseudopotential construction. 
As in the two-particle case we can use the total relative 
angular momentum L of the cluster of TV particles as a 
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useful parameter. Again by rotational invariance of the 
interaction, L will not be changed by the interaction V 
between the particles of the the cluster. However, here 
the parameter L is not sufficient to fully describe the N- 
particle wavefunction as it is in the two body case (See 
Eq. [T|). We must therefore instead write an orthonormal 
set of possible relative wavefunctions \L,q) all having the 
same total relative angular momentum between the N- 
particles. Here the index q runs from 1 to the number of 
states p(N, L) in the basis for that given L and N. We 
will elaborate more on the structure of these wavefunc- 
tions below. 

Given this basis of relative wavefunctions, we define 
TV-particle pseudopotentials by: 

V q L q N = (L,q\V(v 1 ,...,V N )\L,q'). (5) 

As in the two-body case, by translational invariance of 
the interaction, we need not specify the center of mass 
wavefunction in the definition of the pseudopotential. 
Again, by rotational invariance of the interaction, the 
relative angular momentum L must be conserved (the 
matrix clement is diagonal in this variable). However, 
there is no need for the matrix element to be diagonal 
among the states q with the same angular momentum. 
Hence for each TV and L we define a p(N, L)-dimensional 
Hcrmitian pseudopotential matrix with indices q and q'. 

The Hamiltonian can be written in terms of these new 
pseudopotentials via a resolution of the identity (and 
making use of the fact that the interaction potential is 
rotationally invariant by construction): 

H = E E \ l >Q) v l,n{ l >Q'\> ( 6 ) 

z'l<...<ijv L,q,q' 

where the sum over i\, ..An indicates which TV particles 
are considered in a particular term of the sum and it is 
implied that \L,q) is the relative wavefunction for that 
given set of particles. 

We now turn to the issue of determining the basis 
\L,q). These wavefunctions must be made of LLL vari- 
ables, and must be overall symmetric for boson wavefunc- 
tions and antisymmetric for fcrmion wavefunctions. Fur- 
ther, the basis states represent relative motion so they 
should be translationally invariant (i.e., the center of 
mass co-ordinate should not appear). 

To be more specific, let us factor out the ubiqui- 
tous Gaussian exponential factors throughout the dis- 
cussion (indeed, these factors are geometry dependent 
anywa y 10 i 13 ). The remaining wavefunction must be a 
homogenous analytic polynomial of degree L in the vari- 
ables Zi analogous to Eq. [T] Translational invariance of 
the wavefunction implies that the polynomial must be 
invariant under any global shift in all of the co-ordinates 
Zi — > Zi + a for any complex number a. Finally, the wave- 
function must be overall symmetric or antisymmetric for 
bosons or fermions respectively. 

The enumeration of such polynomials is a task that 
turns out to be fairly straightforwar d 13 ' 47 . First we 



note that the problem of enumerating the antisymmetric 
wavefunctions is essentially equivalent to that of enumer- 
ating the symmetric wavefunctions. To see this we note 
that any antisymmetric polynomial in TV variables can be 
written as a symmetric polynomial in TV variables times 
a "Jastrow" factor (Vandermonde determinant): 

N 

ji...n = n (zi- Zj), (7) 

«>j=i;i<j 

Thus there is a precise isomorphism between homoge- 
neous symmetric polynomials of overall degree L and 
homogeneous antisymmetric polynomials of degree L + 
TV(TV- l)/2. 

To establish a complete basis of translationally invari- 
ant symmetric polynomials we use the basis of the ele- 
mentary symmetric polynomials, which are defined in the 
following way: 

e m ,N (zi, z N ) = (8) 

!J2 Zi!—Zi m m<N 

0<i 1 <i 2 <.--<i m <N 
otherwise. 

All possible symmetric polynomials in TV variables can 
be written as sums and products of these generators (i.e., 
these generate the "ring" of symmetric polynomials). 

To impose the condition that the polynomials are 
translationally invariant we shift each variable by the 
overall center of mass coordinate to give 

1 N 

Zi=Zi - — ^Zj. (9) 
j'=i 

Then by writing elementary symmetric polynomials 
of the relative co-ordinates Zi, we obtain generators 
for the ring of translationally invariant symmetric 
polynomia l] 13 ' 47 . It is easy to check that 

ei.iv (zi,z 2 ...zn) = 0. (10) 

therefore there is one fewer generator once we impose 
translational invariance. It can be show n 13 ' 47 that the 
remaining generators e m ^ (zi,Z2 . . . zn) for 1 < m < TV 
do not vanish and are still linearly independent. 

Given that we know the generators, with some combi- 
natorics, we can calculate the dimension d sym (L,N) of 
the space of translationally invariant symmetric polyno- 
mials in TV variables and of degree L. A table of the 
values of c? S ym (L, TV) are listed in Rcf. The ana- 

lytic formula for d sym (L, TV) is reproduced in Eq. ID2I in 
appendix [Dl and these values are identical to the spin 
polarized cases presented in table U (see rows with the 
maximum S'jjQgQjj value, i.e., the top row, for each value 
of TV). 

As an example, consider a translationally invariant 
symmetric polynomial of order L = 4 in N = 4 variables. 
From our allowed generators, the only basis states we can 
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construct of order L = 4 are given by [e 2 ,4 {z\, Z2, Z3, S4)} 2 
and 64,4 (zi,%2, S3, Hence d sym (4, 4) = 2 (c.f tabic [TJ) . 

Our objective was to determine a basis for the states 
\L, q). For bosonic wavefunctions, using combinations 
of these generators, we have found an appropriate basis 
of translationally invariant symmetric polynomial wave- 
functions. The basis states given here are not orthonor- 
mal, but can easily be orthonormalizcd by hand (indeed, 
the concept of orthonormality depends on the integration 
measure. For example, on the sphere it differs from the 
usual Gaussian we are familiar with on the plane). For 
fermionic wavefunctions, we simply multiply these sym- 
metric functions of degree I by a Jastrow factor Ji...jv 
to give a basis for the space of translationally invariant 
antisymmetric polynomials of degree L + (TV(TV — l))/2 
in TV variables. 

Having determined our basis of states \L,q) it is then a 
straightforward matter to construct our pseudopotential 
representation of any given Hamiltonian using Eqs.[5]and 
® 

III. RESULTS FOR MULTI-COMPONENT CASE 

In this section, we present our main results for the 
multi-component case. The essence of our objective in 
this section is, analogously to the single component case, 
to determine how many different linearly independent 
TV-particle wavefunctions might exist with a fixed an- 
gular momentum L, and to form a complete basis for 
these states. For the multicomponcnt case, we may be 
able to classify these states by some additional quantum 
numbers (such as overall spin in the case where we are 
considering the multiple components to be multiple spin 
states). While some amount of formalism is necessary 
in this section, it will be minimized. The more detailed, 
and more formal, derivations are left to the appendices. 

Very generally we will consider an TV-body interaction 
Hamiltonian which may depend on an internal degree of 
freedom Si (such as spin) of each particle 

H= V(v 1 ,8 1 ;...;t n ,s n ). (11) 

il <«2 <■■■<!» 

We will assume that the interaction V is translationally 
invariant, and rotationally invariant in positional space 
(i.e., under rotation of the r variables), but we do not 
necessarily assume the interaction is invariant under any 
particular symmetry of the internal degree of freedom 
(for example, if we are considering particles with spin, 
we do not assume rotational symmetry in spin space). 

Analogously to our approach in the spinless case, to de- 
compose this interaction into pscudopotcntials, we need 
to construct a complete set of states for an TV-particle 
cluster. Again factoring out the center of mass degree 
of freedom of the cluster, let us write a complete set of 
relative wavefunctions as 

\L,q), (12) 



where L is the relative orbital angular momentum of the 
cluster and q enumerates all basis states with this value 
of L (note that the index q indicates not only different 
spatial wavefunctions, but also the different possible con- 
figurations of the internal degree of freedom e.g. spin.) 
Given such a complete basis, we can always define the 
corresponding pscudopotcntials as the matrix elements 
of the form: 

V£$ = (L,q\V(r 1 ,s 1 ;...;r N ,s N )\L,(]'}. (13) 

Note that due to the rotational invariancc of the poten- 
tial, pseudopotentials are always diagonal in L, and due 
to the translational invariance of the potential, the center 
of mass degree of freedom does not appear. This expres- 
sion is analogous to Eq. [5] above. 

We can now use a resolution of the identity to rewrite 
the Hamiltonian, Eq. II H in the following general form: 

H = E E %*)vziM\> ( 14 ) 

ii<...<jjv L.q.q' 

analogous to Eq. EJ 

Thus, our task in this section is simply to determine 
the complete linearly independent basis \L,q) when the 
particles have an internal degree of freedom. 



A. Two Component Case: Spin 1/2, Bilayers, etc. 

We begin with the simplest and most important multi- 
component case — the two-component case. This case 
applies, for example, to spin- 1/2 fermions, such as (un- 
polarizcd) electrons where each fermion has two possible 
internal states (spin- up and spin-down). This case also 
applies to (spin-polarized) bilayer o 28 ' 29 , where the layer 
index (or iso-spin) corresponds to the two-state system. 
We may also consider bosons with two internal states 
which are frequently called "spin- 1/2 bosons" (although 
this nomenclature is not strictly correct). These two 
internal states could be two available hypcrfine states, 
or potentially two possible layers or any other orbital 
indexes. Whatever the origin of the two possible states, 
we will use the nomenclature of spin for simplicity. (Note 
that in the case of more than two components which we 
discuss in section IIII CI below, the language of spins be- 
comes somewhat less useful.) 

In order to describe wavefunctions which depend on a 
spin degree of freedom, it is convenient to work with a 
basis of states which are eigenstates of spin angular mo- 
mentum. These eigenstates are characterized by the spin 
quantum numbers S and S z which are the eigenvalues 
of the combined total spin angular momentum operator 
for TV particles, S 2 , and of the combined z-component 
of spin angular momentum operator for TV particles, S z 
respectively. 

Thus we propose to write a complete basis of states for 
TV spin- 1/2 particles with total orbital relative angular 
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momentum L, total spin angular momentum S 2 and z- 
componcnt of spin angular momentum S z . We denote 
this basis as 

\L,S,S z ,q), (15) 

where q which runs from 1 to the number of states in 
the basis (i.e., the total number of states of N particles 
having L, S 2 and S z ). In the language of Eq. [T21 the 
index q here represents {S, S z ,q}. 

Our first goal will be to determine the dimension of the 
space of the wavefunction basis (the number of q values 
for a given set of spin eigenvalues and a given N, L, S 
and S z ), as has been documented for the spin-polarized 
case (see Eq. |D2] below or Ref. EH). The results of our 
calculation of these parameters for fermionic and bosonic 
cases are shown in table |U For each S, there are always 
25 + 1 different possible values of S z . The table presents 
the number of states for each possible value of S z . (Note 
that some rows of the table are not labeled with a spin 
quantum number Sboson or ■Sfermiorn f° r the bosonic 
or fermionic case, which means that the corresponding 
states cannot occur for a two-component system.) 

Our second goal is to describe the forms of the ba- 
sis wavefunctions. We shall now summarize our results 
leaving the details of the derivation of the forms these 
wavefunctions take to appendices |Bj [C] and [D] below. As 
in the spinless case, the basis wavefunctions, of the form 
\L, S, S z ,q), are equivalent to wavefunctions describing 
a small number of spin- 1/2 particles. Compared to the 
spin-polarized case, these basis wavefunctions are now 
composed of both a spin part and a spatial part. 

Once we have our complete basis, the set of spin- 
pseudopotentials can be defined as in Eq. 1131 

yS,S z , q ;S',S' z , q ' = (16) 

(L, S, S„ q\ V (n, Sl ; . . . ; r N , s N ) \L, S', S' z , q') . 

Using these spin pscudopotentials we can now write down 
an expression for the Hamiltonian Eq. [11] in the general 
form described by Eq. [TU 

As always, rotational invariance in the plane ensures 
that the Hamiltonian is diagonal in L. If the system is 
spin rotationally invariant then the Hamiltonian is also 
diagonal in the eigenvalues S and S z . For a more gen- 
eral interaction however, the Hamiltonian might not be 
spin-rotationally invariant, but nonetheless it is still con- 
venient to decompose the interaction using the spin basis. 

We now turn to the explicit construction of our com- 
plete set \L, S, S z ,q). A wavefunction for fermions must 
be overall antisymmetric whereas a wavefunction for 
bosons must be overall symmetric. However, the wave- 
functions we consider are a combination of both a spin 
and a spatial part, and only the combination of the two 
parts needs to have the overall fermionic or bosonic sym- 
metry. The spatial and spin parts of the wavefunction 
can have more complicated symmetry as long as the two 
parts are appropriately sewn together and the combina- 
tion has the correct overall symmetry. In fact there is 



a direct correspondence between the type of symmetry 
and the spin quantum number S. The more mathemat- 
ical structure of forming this combination is discussed 
in detail in appendix [B] below. Here, however, we will 
simply present the results of this procedure. 

The presentation of our results is divided into two 
parts: firstly we shall introduce primitive polynomials, 
which are the lowest degree polynomials corresponding 
to a particular type of spatial symmetry type for a set 
of TV particles; secondly we shall describe how to use 
these primitive polynomials to construct a spatial func- 
tion of arbitrary degree that still corresponds to a partic- 
ular symmetry type. Each allowed spatial symmetry cor- 
responds to a particular spin eigenvalue for the N parti- 
cles. Merging a spatial wavefunction of a given symmetry 
type (a given spin quantum number) with a correspond- 
ing spin wavefunction will give an overall wavefunction 
with the appropriate fermionic or bosonic symmetry In 
a moment we will discuss how the spatial and spin wave- 
functions are merged. 

In tables ITTI and IIIII we list what we have termed prim- 
itive polynomials, /3l, of degree L. These polynomials 
are written in terms of symmetric polynomials of the 
form of Eq. [5] in terms of the relative co-ordinates as 
in Eq. |H1 Now, however we frequently need to write sym- 
metric polynomials in fewer than all N of the z variables. 
For compactness we shall use the short-hand notation: 

e m,?iZ2^3...^ p = &m,p i , , ...Zip) . (1^) 

Note that we will always take z,; to have the center of mass 
of all N particles subtracted off (as in Eq. [9J independent 
of the value of p. So for example 

e2,245 = e 2 ,3(^2, 24, h) = + z 2 h + Zih- 

We also use the short-hand notation for Jastow factors 
Jn.., P = II %-Zi), (18) 

k . . .i p } ;k<l 

so that for example 

Ji34 = (Si - z 3 )(Si - z 4 ){z 3 - 2 4 ). 

Multiplying an A-particle spatial wavefunction by any 
fully symmetric translationally invariant polynomial docs 
not change its symmetry type. One can thus construct 
spatial wavefunctions with a given symmetry type and 
relative angular momentum L by multiplying a primitive 
polynomial by any fully symmetric (in all N variables) 
translationally invariant polynomial such that the com- 
bined polynomial order of the resultant product is L. 
Thus the most general spatial wavefunction of a partic- 
ular symmetry and a particular relative angular momen- 
tum L is a linear combination of these products of trans- 
lationally invariant polynomials times primitive polyno- 
mials giving a homogeneous polynomial of degree L. 

We shall now demonstrate how to use the primitive 
polynomials via a simple example. Consider a 3-particle 
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electron wavefunction with total spin eigenvalue S = 1/2. 
For relative angular momentum eigenvalue L = 1 we con- 
struct from the relevant entry in table |TT]the only possible 
order 1 polynomial conforming to this type of symmetry, 
namely 

fh = J 12 = {z x - z 2 ) . (19) 

For L = 2 we can use the information in table [TT] to 
construct the valid order 2 polynomials conforming to 
this type of symmetry. Note that due to equation (fTU)) 
61,123 = 0, and so we cannot construct an order 2 polyno- 
mial from J12 (which is order 1) multiplied by any order 
1 translationally invariant fully symmetric polynomial. 
The only possibility according to table [TT] is the second 
primitive polynomial: 

fh = J\i&\si = J12 {zi + h) ■ 

For L — 3 we can have only the polynomial given by 
multiplying the L = 1 result by a fully symmetric trans- 
lationally invariant polynomial of order 2 

(zi - z 2 ) e 2: i23- 

For L — 4 we can have either the L = 1 primitive 
polynomial multiplied by a translationally invariant fully 
symmetric polynomial of degree 3 or the L = 2 primitive 
polynomial multiplied by a translationally invariant fully 
symmetric polynomial of degree 2. The most general 
result is a linear combination: 

(zi - z 2 ) (^4ie 3i i23 + A 2 (zi + z 2 ) e 2 ,i 23 ) . 

or equivalently we have two linearly independent basis 
states in the space of spatial wavefunctions. At each 
polynomial order the number of linearly independent ba- 
sis vectors appearing in these polynomials is precisely the 
dimension appearing in tablclU in this case 0,1,1,1,2... for 
polynomial orders L = 0,1,2,3,4... and so on. Since all 
fully symmetric translationally invariant polynomials are 
generated by e 2; i2 3 and 63,123, the most general spatial 
wavefunction with relative angular momentum L is thus 

(Zl ~ Z2) I E ^/m4,123 e ™123 + 



(Zl+Z 2 ) E ^'m'4',123 e ™123 \ > 
2Z'+3m' 
= L-2 

where the total number of coefficients Ai m and A\, m , ap- 
pearing in the wavefunction takes the value given in table 
III 

To generate the complete basis wavefunctions we must 
combine the spatial wavefunction with a spin wavefunc- 
tion. We shall define a 'primitive spin wavefunction $j of 
a many particle system to be an eigenfunction of the S z 
operator of every particle in the system. For example, 



Ititti) is a primitive spin wavefunction for 5 particles. 
We imagine ordering all of the primitive spin wavefunc- 
tions in lexicographical order (an alphabetical, or last 
letter sequence ordering scheme), so that the first primi- 
tive spin wavefunction (first in the sense of this ordering 
scheme) is given by 

^(N,s z ) = \m -m), (20) 

where the number of spin-up particles minus the number 
of spin-down particles times 1/2 is equal to the total S z 
eigenvalue: 

±(N t -N i ) = S z . 

To combine our spatial wavefunction with a spin wave- 
function, we take a spatial wavefunction corresponding to 
the spin eigenvalue S (i.e., a primitive polynomial in ta- 
ble |TT] or IIIII times any translationally invariant symmet- 
ric polynomial), then multiply this by the first primitive 
spin wavefunction with any valid S z eigenvalue for that 
particular S eigenvalue. To recover the complete basis 
wavefunction we simply antisymmetrize (for fermions) or 
symmetrize (for bosons) the resulting combination in all 
of the co-ordinates. 

Continuing with our example let us consider the case 
of N = 3 where the spin eigenvalue is S z = 1/2. At or- 
der 1, for example, we multiply the spatial wavefunction 
Eq. [TO] with the first primitive spin wavefunction |ttl) to 
construct 

im) (zi - h) ■ (21) 

Then by antisymmetrising in all the particle co-ordinates 
we find the three body wavefunction: 

\L=1,S= 1/2, S z = 1/2) = 2 |m> (zi - h) (00 , 

-2 mt) (zi - h) + 2 iitt) (h - h) ■ [ ' 

Note that we do not write an additional q index here 
(as in Eq. I15p since there is a unique wavefunction with 
this value of N, L, S, S z . Note also that this wavefunc- 
tion is not normalized (as normalization depends on an 
inner product which depends on whether the geometry 
is on the sphere, plane, or torus ). We can obtain a sim- 
ilar wavefunction with S z — —1/2 by simply following 
the same procedure with the primitive spin wavefunction 
An identical result is obtained by instead applying 
the spin lowering operator S~ = E; <ST ^° Eq- E3 

To give another example briefly, to create a bosonic 
wavefunction of four particles with angular momentum 
S = S z = 1 we look on table|nl]under N = 4, S = 1. The 
listed spatial polynomials (such as ei.123 = Zi + z 2 + 03) 
can be multiplied by any overall symmetric polynomial 
(which does not change their symmetry). The spatial 
wavefunction must be then multiplied by the appropriate 
first primitive spin wavefunction |tttl), an( l finally the 
entire wavefunction should be fully symmetrized over all 
particles. 
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The prescription laid out here is sufficient to construct 
a complete basis of states of a two component system 
of TV-particles having overall relative (orbital) angular 
momentum L and spin quantum numbers S and S z . For 
many readers this is the most important result of the 
current paper. 

B. Interlude: Primer on Symmetry Types 

In the above section we explained briefly how we take 
spatial and spin wavefunctions of certain symmetry types 
and sew them together to make a wavefunction with over- 
all bosonic or fermionic symmetry. In this section we 
briefly explain in more detail precisely what we mean by 
a symmetry type. While it is not so essential to under- 
stand this material in order to follow the construction 
of the previous section, in the more complicated cases 
of the next two sections where there are more than two 
components to the wavefunction (spin 1, 3/2, graphene 
with spin and valley, etc) it becomes more crucial to un- 
derstand the concept of symmetry type. The informa- 
tion contained in this section is standard mathematics, 
reviewed here briefly for convenience. See, for example, 
Ref. |H for a more detailed discussion. 

A key intellection which we shall use to describe our 
basis wavefunctions is the concept of a symmetry type 
or antisymmetry type of a wavefunction (the original 
method was introduced in Ref. l49h . These symmetry 
types are intimately related to irreducible representa- 
tions of the symmetric group (the group of permuta- 
tions of N objects). Such irreducible representations 
are visualised by Young frames (closely related to Young 
tableaux) , a series of N boxes arranged in left -justified 
rows and columns (a more detailed description of the rep- 
resentation theory of the symmetric group and of Young 
tableaux is provided in appendix [2]) • We shall denote 
the shape of a frame containing Ni boxes in the ith row 
by A = [Ni,N 2 ,...] with Ni > Nj > if i < j and 
TV = J^i Ni ■ Thus the shape of each frame represents a 
unique partition of the integer N into pieces [N\, N2, ■ ■ ■]■ 

We will use the term conjugate to refer to Young 
frames related by a reflection along the diagonal, for ex- 
ample, the following frame of shape A = [3, 1]: 



is conjugate to the following frame of shape A = [2, 1, 1 



Representations of the symmetric group associated with 
conjugate Young frames are called contragradient rep- 
resentations. We will use the notation A to denote the 
conjugate Young frames (or integer partition) to A. 

We now explain how Young frames can be used to de- 
scribe the permutation symmetry of a function. Given 



a Young frame A = [N\, iVa, • • ■], being an integer par- 
tition of N, we say a function has normal form with 
symmetry A if (a) it is symmetric in the first N\ vari- 
ables, symmetric in the next N2 variables, symmetric in 
the next N3 variables and so forth and (b) we cannot 
use permutation and linear combinations to construct a 
nonvanishing function which is symmetric in greater than 
Ni variables, we cannot use permutation of the variables 
Ni + 1, . . . N and linear combination to construct a non- 
vanishing function which is symmetric in greater than 
N2 of these variables, we cannot use permutation of the 
variables Ni + N2 + 1, . . . N to construct a nonvanish- 
ing function symmetric in greater than N3 of these vari- 
ables and so forth. In other words, permutations and 
linear combinations cannot make the function any more 
symmetric than it already is. By permutation and linear 
combination any function can be brought to normal form 
to clearly show its symmetry type (the detailed method 
of reduction to normal form is given in Ref. l48h . 

Antisymmetry type is entirely analogous except that 
one should replace the word symmetric with antisym- 
metric in the above definition. A nontrivial theorem^ 
is that any function which has symmetry type A has the 
conjugate antisymmetry type, denoted by A. 

One can check that the polynomials given in Tables 
im and |TTT] are listed by symmetry type. For example, 
in table IIII1 the N = 5 particle functions listed of type 
[3, 2] are clearly symmetric in the first three variables 
(zi,Z2,Z3) and then also symmetric in last second two 
variables (2:4,2:5). Further, any attempt to symmetrize 
in four variables will vanish. Thus these entries are in 
normal form. For table ITTI the symmetry type is less ob- 
vious (it is not in symmetric normal form), but the an- 
tisymmetry type is evident (the polynomials are listed 
in antisymmetric normal form). For example, symmetry 
type A = [2,2,1] must be equivalent to antisymmetry 
type [3,2]. One can check that the polynomials labeled 
by [2,2,1] in table [TT| are indeed antisymmetric in the 
first three variables and then also antisymmetric in the 
last two variables and attempts to antisymmetrize in four 
variables will fail. 

Having defined symmetry type, the general procedure 
to create a wavefunction with bosonic symmetry involves 
sewing together a spin and a spatial wavefunction of the 
same symmetry type. Similarly to create wavefunctions 
of fermionic symmetry one must sew together spin and 
spatial wavefunctions with conjugate symmetries. The 
details of this procedure are presented in appendix [5] be- 
low. Here we shall simply state a few key results derived 
in that appendix. The first result is that, for a two- 
component system, the spin wavefunctions of a particu- 
lar symmetry type are in one-to-one correspondence with 
spin eigenfunctions, that is, such functions are eigen- 
functions of the 5 2 operator. For example the follow- 
ing function is an cigenfunction of S 2 with eigenvalue 
S = S Z = 1/2: 

4. (2 im> - mt> - iitt» • 
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This can be shown by explicitly applying an S 2 oper- 
ator. One can also check that this function has sym- 
metry type [2, 1], since it is symmetric in the first two 
arguments. In order to generate a fully antisymmetric 
combined spin and spatial wavefunction, the procedure 
is to combine the above spin wavefunction with a spa- 
tial wavefunction of conjugate symmetry type, and then 
antisymmerise the resulting construction. To generate 
a symmetric combined function, we associate that spin 
cigcnfunction with a spatial function of the same sym- 
metry type and then symmetrize that construction. If, 
for example, we wanted to generate a 3-electron wave- 
function then the corresponding spatial part would be a 
polynomial of conjugate symmetry type [2, 1] (note that 
[2, 1] is self-conjugate). 

The second key result shown in appendix [C] is that, 
for a two-component system, the result of the construc- 
tion procedure described in the previous paragraph can 
be obtained by using just the first primitive spin wave- 
function (defined in Eq. I2"0"f instead of the full spin 
cigcnfunction. The simplest construction procedure is 
then to associate a spatial function of a given symme- 
try type with a first primitive spin wavefunction only, 
and then to (anti) symmetrize or symmetrize that result 
to obtain an overall (anti)symmetric combined spin and 
spatial wavefunction. The reason for this simplification 
lies with the (anti)symmetrisation operation: as long as 
the spatial part corresponds to a particular symmetry 
type, the the (anti)symmctrisation procedure will auto- 
matically impose the correct symmetry type on the spin 
part. 

Once we have more than two-components in the wave- 
function (higher spin, or multiple spins and valley de- 
grees of freedom) this simplification no longer fully holds, 
nevertheless our guiding principle is to try to specify a 
(generalized) spin wavefunction of the simplest possible 
form, which can be associated with a spatial wavefunc- 
tion of a particular symmetry such that the result of 
(anti)symmctrising that construction is an appropriate 
basis wavefunction. 

A system built from N particles with n multiple com- 
ponents is described mathematically by the irreducible 
representations of the Lie algebra of the group SU(n). 
For example, we refer to spin- 1/2 particles as being de- 
scribed by representations of the Lie algebra of SU(2). Ir- 
reducible representations of SU(n) are in correspondence 
with a subset of irreducible representations of the sym- 
metric group of N objects^ and hence the symmetry 
types. The simplest way to visualize this correspondence 
is via Young frames associated with these irreducible rep- 
resentations. The set of Young frames describing the irre- 
ducible representations of SU(n) corresponds precisely to 
the subset of symmetric group Young frames containing 
no more than n rows; the corresponding set of conjugate 
Young frames contain no more than n columns. There 
is a corresponding restriction on the possible symmetry 
types which can be used, for example in the case of a two 
component system the possible symmetries of the spin 



and spatial wavefunction are restricted to being symmet- 
ric or antisymmetric in two subsets of particle indices. 
Hence in table IIIII where bosons are considered the only 
allowed symmetry type for the two-component case (spin- 
1/2 case) are those where there are only two columns 
of the frame A, whereas in table [TT| where fermions are 
considered the only allowed symmetry type for the two- 
component case (spin- 1/2 case) are those where there are 
only two rows of the frame A (or two columns of the con- 
jugate frames). 

C. Extension to Systems with n Components 

Having described the general properties of general 
symmetry types, we shall now turn our attention to sys- 
tems containing in general n components. In this case the 
possible symmetries of the spin and spatial wavefunctions 
are symmetric or antisymmetric in up to n subsets of 
particle indices. The corresponding symmetry types are 
A = [Ni,N 2 , ■ ■■,N n ] (for bosons) or conjugate symmetry 
types A (for fermions). In order to generate a complete 
basis for interactions with multi-component symmetry, 
we have chosen to decompose the inter-particle inter- 
action into a basis of all the possible combinations of 
symmetry types that the spatial and spin wavefunctions 
can have. The forms of the spatial wavefunctions remain 
translationally invariant analytic polynomials. 

In table U we have listed the dimensions of the space of 
polynomials for an n-component system (n possible val- 
ues of the internal degree of freedom) for TV < 5 particles. 
The dimensions are labeled by the symmetry type A of 
the corresponding pseudopotential basis functions. Ad- 
ditionally, primitive polynomials for all symmetry types 
with N < 5 are listed in tables |H] and HID We stress 
that not all symmetry types are available to a given n- 
componcnt system — such a system can only be classified 
into symmetry types of the form A = [Ni, N2, N n ] (for 
bosons) or conjugate symmetry types A (for fermions). 
In particular, for the two-component case, one can see 
that only certain symmetry types are allowed depending 
on whether the system is bosonic or fermionic (e.g., the 
fully symmetric [3] type is not allowed for three fermions). 

In order to construct complete basis wavefunctions, we 
also require a modified "spin cigcnfunction" to describe 
the multiple internal degrees of freedom. We shall refer 
to this wavefunction as a generalized spin wavefunction. 
When a particle has an internal degree of freedom s, let 
us notate the different possible values of this degree of 
freedom by a, j3, 7 and so on. The notation can be 
interpreted, for example, in the two-component case of 
spins as, a = \f) and j3 = \\); but now there n possible 
values that our parameter can take. For a system with 
multiple internal degrees of freedom, we similarly define 
primitive generalized spin wavefunctions such as 

d = \a(P/3a) . (23) 

It is useful to define the index A 2 which represents the 
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list {N a! Np, ...,N^} of how many times each type each 
component {a., ft, . . .) occurs in the primitive generalized 
spin wavefunction. So for example, in the primitive gen- 
eralized spin wavefunction 1 0:0:7) for three particles we 
have the index X z = {2, 0, 1}. 

A consequence of the increased number of component 
degrees of freedom is that the complete basis wavefunc- 
tion can no longer be generated solely by using a so-called 
"first" primitive generalized spin wavefunction in analogy 
with the two-component case. Instead, we consult a list 
of the primitive components that occur for each sym- 
metry type in order to generate a complete basis. These 
results are presented for the case of three-component and 
four-component systems in table IIV1 

As an example, consider the three-component N = 3, 
A = [2,1] and X z = {1,1,1} state. From table [TV] we 
see that there is a choice of two possible primitive gen- 
eralized spin wavefunctions, \a/3j) or \aj(3). Using these 
functions to construct a complete basis wavefunction will 
generate two orthogonal wavefunctions. (An interesting 
aside is that the three-component, N = 3 generalized 
spin wavefunctions arc identical to the color wavefunc- 
tions describing baryons in the SU(3) quark model, see 
for example Ref. IHoh 

The problem of constructing n-component iV-particlc 
wavefunctions with prescribed symmetry is equivalent to 
the problem of constructing tensor products of N funda- 
mental irreducible representations of the Lie algebra of 
SU(n). We shall explain the mathematical details of that 
procedure in appendix [A] Our general aim is to decom- 
pose our parameter space of wavefunctions into a basis 
of wavefunctions labeled by irreducible representations of 
the Lie algebra of SU(n). We shall describe the general 
method for selecting the appropriate primitive general- 
ized spin wavefunctions in appendix [C] Due to symme- 
try arguments, note that it is not possible to construct 
linear combinations of certain types of primitive general- 
ized spin wavefunction corresponding to certain symme- 
try types. For example, it is clearly not possible to gener- 
ate a fully antisymmetric generalized spin wavefunction 
from |o...o). 

The method to construct a full basis wavefunction is 
similar to the two-component case: first we choose a 
primitive generalized spin wavefunction from table IIVI 
corresponding to a given symmetry type A and having a 
given X z . For a boson wavefunction we construct a poly- 
nomial corresponding to a symmetry A and at a given 
order L using the primitive polynomials in table Hill The 
full basis wavefunction is generated by putting these two 
parts together and then fully symmetrizing. To obtain a 
fermion wavefunction we would instead construct a poly- 
nomial corresponding to a symmetry type A using the 
primitive polynomials given in table UH and the complete 
basis function is given by fully antisymmctrizing. We 
shall explain how these results were derived in appen- 
dices IB O and O 

To obtain the total size p(N, L) of the space of pseu- 
dopotcntials, \L, q), for a given L, we need to include both 



the number of polynomials p(N, L, A) with a given sym- 
metry type A and the number of generalized spin wave- 
functions for that same A but with different values of X z 
(for example, in the two-component case the spin wave- 
functions have 25+1 values of S z ). For clarity, let us dis- 
cuss an example: a system of TV = 4 bosons with n = 3 
components. Say that we want to construct the basis 
space of wavefunctions for an interaction at L = 4. The 
general form of the basis wavefunctions is \L = 4, q). In 
this case, the possible symmetry types that the wavefunc- 
tion can have are A = [4], [3, 1], [2, 2] and [2, 1, 1]. Look- 
ing at table HI the wavefunctions with these symmetry 
types have the respective polynomial space dimensions 
of p(4,4, [4]) = 2, p(4,4, [3, 1]) = 2, p(4,4, [2,2]) = 2 and 
p(4, 4, [2, 1, 1]) = 1. Further, the corresponding dimen- 
sions of the space of generalized spin wavefunctions are 
given by counting all of the primitive generalized spin 
wavefunctions in table ITVl which arc respectively 15, 15, 
6 and 3. The total dimension of the space of L — 4 pseu- 
dopotentials, that is, the total number, p(N, L), of possi- 
ble values of q, is thus 15x2 + 15x2 + 6x2 + 3x1 = 75. 
Each of these basis wavefunctions is obtained by multi- 
plying the spatial part of the wavefunction by the gener- 
alized spin wavefunction and symmetrizing (for bosons in 
this case). The pseudopotential matrix at L = 4 would 
then be a 75 dimensional Hcrmitian matrix! 

In an system with n-components, if the Hamilto- 
nian has full SU(n) symmetry, then the pseudopotential 
Hamiltonian, Eq. 1141 will be diagonal in both A and X z 
variables. 



D. Systems with Spin Rotational Symmetry 

Often when we consider particles with n internal states 
(components), these multiple states may have arisen from 
particles that had spin J — (n — l)/2. If the Hamilto- 
nian is spin-rotationally invariant, then it is very useful 
to consider a set of basis wavefunctions that are eigen- 
functions of S 2 and S z such that the Hamiltonian Eq. 03] 
is diagonal in these variables. In fact, even if the Hamil- 
tonian is not quite spin-rotationally invariant, it may be 
very convenient to work in such a spin basis. We thus 
would like to construct basis states which correspond to 
both a given symmetry type and to a given S 2 eigenvalue 
via linear combinations of the primitive generalized spin 
wavefunctions given in table IIVI Mathematically speak- 
ing, this method corresponds to the decomposition of rep- 
resentations of the group SU(n) into representations of 
SU(2); the procedure for doing the decomposition is de- 
scribed in Ref. |48|. For a system constructed from spins 
of magnitude J, it is to be understood that our notation 
should be interpreted as a corresponding to the highest 
spin state J, (3 corresponding to spin J— 1 and so on. The 
total z-component of spin angular momentum eigenvalue 
is given by: 

S z = JN a + (J - l)Np + (J - 2)/V 7 ... - JN ( 



12 



For spin-1 systems, appropriate linear combinations of 
the primitive generalized spin wavefvmctions forming 
spin-1 eigenfunctions are listed in table [Vj Note that 
these spin eigenfunctions are listed such that S = S z , 
and so additional spin eigenfunctions can be constructed 
by explicitly lowering the S z eigenvalues using the S~ 
operator. Due to the increased spin degree of freedom, 
we find that there can now be several values of S asso- 
ciated with each symmetry type A (by contrast in the 
spin-1/2 case where the spin eigenvalue S corresponds to 
a unique symmetry type). In general it is also possible 
for multiple occurrences of the same S value to occur for 
a given symmetry type, however this feature only occurs 
for N > 5. For systems constructed from spin-3/2 par- 
ticles the procedure is along similar lines and we shall 
outline the general method for constructing such spin 
eigenfunctions for arbitrary spin in appendix [C] 

The method to construct a basis wavefunction is once 
again the same as in the two-component case. This time 
we start by choosing a one of the linear combinations 
of primitive generalized spin wavefunctions from table 
[V] corresponding to a given A and a given S 2 eigenvalue. 
The linear combination of generalized spin wavefunctions 
will be associated with a polynomial corresponding to 
that symmetry A for bosons (or conjugate symmetry A 
in the case of fermions) and at a given order L. The 
complete basis wavefunction is given by symmetrizing (or 
antisymmetrizing) the product of spin and spatial parts. 
We shall explain how these results were derived in ap- 
pendices [B] and [C] 

E. Spin and Valley: the Case of Graphene 

More generally, a particle may have several different 
internal degrees of freedom. A particularly important 
example of this is graphene which has two spin states and 
two valley states resulting in four possible internal states. 
The Hamiltonian may have full SU(4) symmetry^— or, 
for example, it could be symmetric only under the SU(2) 
rotations of the spin. 

It is useful always to think of this in terms of a full 
SU(4) symmetry (which may be broken). We have al- 
ready enumerated all possible spatial wavefunctions with 
symmetry types corresponding to SU(4) (for N < 5): the 
dimensions of the corresponding polynomial spaces are all 
contained within table HI and corresponding spatial wave- 
functions for these symmetries are in tables [TT] and |TTT] re- 
spectively for bosons and fermions respectively However, 
in the case where the system has an SU(2) spin rotational 
symmetry in each valley (but not the full SU(4) symme- 
try) it is useful to decompose these states into their S 2 
eigenstates so as to exploit the symmetries of the prob- 
lem as much as possible. Mathematically speaking we are 



decomposing SU(4) representations into representations 
of SU(2) x SU(2). 

In our generalized spin wavefunction notation, it is to 
be understood that in this situation we interpret a as 
corresponding to spin-up in the first valley, f3 as corre- 
sponding to spin-down in the first valley, and similarly 
7 and 5 as representing spin-up and spin-down states in 
the second valley. The total z-component of spin angular 
momentum eigenvalue is thus given by: 

S, = l{N a + Nj-±(N fi +N s ). 

We can then construct spin eigenfunctions by choosing 
appropriate linear combinations of primitive generalized 
spin wavefunctions. In this case there are necessarily 
multiple linearly independent S 2 eigenstates correspond- 
ing to each SU(4) symmetry type. The appropriate linear 
combinations of primitive generalized spin wavefunctions 
are given in table IVD 

To clarify: each primitive generalized spin wavefunc- 
tion is labeled by A, S, and S z , and, as shown in the table, 
there may be several basis wavefunctions with the same 
value of all of these variables and we can distinguish be- 
tween these cases by introducing a further index k. Note 
that in the table only states with S — S z arc shown; 
other values of S z are obtained by applying a lowering 
operator. 

Given any primitive generalized spin wavefunction in 
table IVII we follow the same procedure as in the above 
cases — for bosons (fermions) we multiply it this primi- 
tive generalized spin wavefunction by a polynomial with 
the same (conjugate) symmetry type A (A) and then 
(anti)symmctrize over all variables. 

In appendix [C] we shall explain how we derived the 
spin eigenfunctions in table IVTl In that appendix we shall 
also explain further possible extensions to the cases given 
here. For example, we shall explain how to go about 
constructing an appropriate pseudopotcntial basis to de- 
scribe systems where there are more than two interacting 
degrees of freedom per particle, or systems where the de- 
grees of freedom themselves have more than two possible 
values (i.e., are higher spin than 1/2). 



F. Tables 

All of the results presented in the tables here have 
been generated using a Mathematica™-based computer 
program^. The program is capable of producing spatial 
and generalized spin wavefunctions for systems of up to 
N = 5 and for an arbitrary number of components, al- 
though we only have space to give a selection of results 
for the most relevant cases here. 
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TABLE I: Dimensions of the polynomial spaces: the number of independent parameters occurring in the spatial part 
of the wavefunction at order L for quantum Hall states containing TV particles and with symmetry type A. Where 

applicable to the two-component case, which is described in terms of spin cigenfunctions, we list the spin eigenvalues 

S corresponding to either a fermionic (Sfermion) or a bosonic (S'boson) wavefunction. For each S listed there will be 
25+1 possible values of S z . The symmetry types labeled by A are introduced in section UlI Bl We explain in the 

text that not all symmetry types A correspond to a spin eigenvalue Sf erm j on or Sb oson , and we indicate these cases 

by a - in the table. 
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TABLE II: Primitive polynomials /3l of degree L for construction of fermion wavefunctions, listed by their symmetry 
type A. Note that for compactness we use the notation e m> i l i 2 ...i defined in Eq. [T7|to be the elementary symmetric 
polynomial of degree m in the p variables z^, . . . zi , and we similarly write Jastrow factors on p variables as i 
as defined in Eq. 1181 The column labeled Sfermion ^ s m reference to the case of spin-1/2 (two component) particles 
only — entries labeled with a - do not occur for the spin-1/2 case. The symmetry type A pertains more generally. 



(a) N 



X 


'-'fermion 


Primitive Polynomials 


[2] 





1 


[1,1] 


1 


Jl2 



(b)N 



X 


'-'fermion 


Primitive Polynomials 


[3] 




1 


[2,1] 


1/2 


Pi = Jia 

P2 = Jia(ei,ia) 


[1,1,1] 


3/2 


J VIZ 



(c) JV 



A 


'-'fermion 


Primitive Polynomials 


[4] 




1 


[3,1] 




Pi = Jl2 

P2 = Jia(ei,ia) 

Pz = J12 (35 3 ei,i2 + e?, 12 - e 2 ,i2 + 3z|) 


[2,2] 





Pi = J12J34 

Pa = Jia^34(ei l i2) 


[2, 1, 1] 


1 


Pz = Jl23 

Pa = ^123(61,123) 
Pb = (62,123) 


[1,1,1,1] 


2 


J1234 



(d) N = 5 



A 




fermion 


Primitive Polynomials 


[5] 




1 


[4, 1] 




Pi = J12 

P2 = Ji2(ei,ia) 

P3 = J12 (3ef >12 - 3e 2 ,i2 + 2e 2 ,34s) 

Pa = J12 (5ei 12 — 4e 2 ,i2ei,i 2 — 6e2,34sei,i2 - 6e3, 3 4s) 


[3,2] 




P2 = J12J34 

Pz = J12 J34 (ei,i2 + £1,34) 

Pa = J12 J34 (ef,i2 + 2ei,34ei,i2 + ef i34 — e2,i2 — e 2 , 3 4) 

Ph = J12J34 (ei,i2 + 5ei i 34ei jl2 + 5ei, 34 ei,i2 + 4e 2 ,i2ei,i2 + 62,3461,12 + ei, 34 + ei, 3462, 12 + 4ei,34e2,34) 
/ 261,12 + 661,3461,12 + 8ei i34 ei. 12 — e2.i2ei. 12 — 2e2,34ei : i2 + 6e 134 ei, 12 + 2e 2 ,i2e2, 34 \ 

P6 = J12J34 4 2 2 2 2 

V — ei, 3 4e2,i 2 ei,i2 — ei, 3462,3461, 12 + 2ei, 34 + 2e 2 ,i2 + 2e 2 ,3 4 — 2ei i 3 4 e2,i 2 — 61,3462,34 / 
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TABLE II: Continued 
(e) N = 5 (cont.) 



A 


"^fermion 


Primitive Polynomials 


[3, 1, 1] 




Pi = Jl23 

04 = Jl23 (ei,123) 

#3 = Jl23 (-Aie? )123 - 2,4162,45 + 2yl 2 e2,123) 

A> = ^123 (3e5,i23 + 5e 2 ,45ei,i 2 3 — e2, 12361,123 — 263,123) 

Pt = J123 (5ef, 12 3 — 6e2,45ef, 12 3 — 3e3, 12361, 123 + 6e|, 45 — 662,4562,123) 


[2,2,1] 


1/2 


Pi = J123J45 

Pb = >7l23 ^45 (ei,123) 

/3e = J123J45 (ei,i23 + e2,45 — 62,123) 

Pr = Jl23-/45 (3ei,i 2 3 — 2e2,45ei,123 + 63,123) 

Pa = J123 J45 (et,i23 + 3e2,45ei,i23 + 62,12361,123 — 63,12361,123 + e 2 ,i 23 — 62,4562,123) 


[2,1,1,1] 


3/2 


Pd = Jl234 
p7 = Jl234 (ei,l234) 
Ps = Jl234 (ei,i234) 
P9 = Jl234 (ef,i234) 


[1,1,1,1,1] 


5/2 


J12345 
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TABLE III: Primitive polynomials /3l of degree L for construction of boson wavefunctions, listed by their symmetry 
type A. As in Table UTT the symbol e m ,iii 2 ---i p is defined in Eq. ITTland Ji 1 ....i v is defined in Eq. [T5] The column labeled 
£>boson is in reference to the case of spin-1/2 (two component) particles only — entries labeled with a - do not occur 
for the spin-1/2 case. The symmetry type A pertains more generally. 



(a) AT 



(b) N 



A 


"^boson 


Primitive Polynomials 


[2] 


1 


1 


[1,1] 





Jl2 



A 


'-'boson 


Primitive Polynomials 


[3] 


3/2 


1 


[2,1] 


1/2 


01 = ei,i2 

02 = e?,i2 + 2e 2 ,i2 


[1,1,1] 




Jl23 



(c) N 



A 


'-'boson 


Primitive Polynomials 


[4] 


2 


1 


[3,1] 


1 


01 = ei,i23 

02 = (e?,123 + e2,123) 

03 = (ei,123e2,123 + 3e3,i23) 


[2,2] 





Pi = (e?,i2 + 2e 2 ,i2 + 2e 2 , 3 4) 

Pi = (6^12 + 4e 2i i2e? j:L 2 + 462,346^12 — 261,12 — 2el,3 4 — 12e2, 1262,34) 


[2, 1, 1] 




Pz = 32|ei,i2 - 2361,12 + 22 3 e2,i2 - ef,i2 + ei,i2e 2 ,i2 + 223 
Pi = 2361,12 — 32|e?,i2 — 2 3 ef,i2 + 62361,1262,12 + 3e?,i 2 e2,i2 

05 = 2 3 eli2 + 32|e?,i2 + 22|ef,i2 - 42 3 e2,i2ei,i 2 - 62362,1261,12 - 42 3 e|,i 2 - 42|e 2 ,i2 - e 2 ,i2ei,i 2 - 2e!,i 2 ei,i 2 


[1,1,1,1] 




J1234 



(d) N = 5 



'boson 



Primitive Polynomials 



[5] 



5/2 



Pi = ei,1234 



[4, 1] 



3/2 



P2 = 

03 = 

04 = 



3ef,i2 3 4 + 2e2,1234) 

5ei,i2 34 — 3e2,1234ei,l234 + 3e3,i234) 

3ei,i234 + 4e2,1234ei,i234 — 4e 3 ,l234ei,l234 — 2e2,i234 + 4e4,i234) 



[3,2] 



1/2 



P2 = 

03 = 

04 = 

05 = 

06 = 



2e3,i2 3 ei,i23 — 3el 45 — 123 



12e2, 4562,123) 



e?,123 + 3e2,45 + £2,123) 
ei, 12362, 45 + ei, 12362, 123 + 3e3,i2 3 ) 

e l,123 + 6e2,45e?,i23 + 2e2,123e?,i23 
e2,45ei,l23 + e2,123ei,123 + 7e3,l23ei,l23 

^ 2ef,i23 + 18e 2 ,45ef,i23 + 6e 2 , 12361,123 — 6 e 3 ,i23ef,i 2 3 — 27el, 45 e?, 12 3 - 9e|, 123 ei,i 2 3 ^ 

-60e2,45e2,i23ei,i23 — 60e2, 4563,123 ei,i23 + 12e2, 12363, 12361,123 + 6el, 4 5 + 2 e|,i23 
y +60e2,45e|,i23 - 3e3,i23 + 60e|, 45 e 2 ,i23 



e 2,45 e l,123 — e2,123 e l,123 — 12e2,45C3,123 — 3e2, 12363,123) 



Pz = J45 (2ei,i23 + 3e2,45 — 2e2,i23) 

Pi = J45 (2ei,i23e2,45 + 3ei,i23e2,123 — 3e3,123) 

^ 20ylie2,45e?,i23 - 2A 2 e2, 4 5ef,i 2 3 + 20ylie2,i23ef ,123 + 2yl 2 e 2 , 1236?, 123 - 4yl 5 e 3 , 12361,123 

06 = ^45 — 10Aie2,4 5 ei,i23 — 3v42el, 45 ei,i23 — 20Aiel, 12 3ei,i23 — 2A2e|, 12361, 123 — 4(L4ie2, 4562, 12361,123 

y +4yl2e2, 4562, 12361, 123 + 40Aie2, 12363, 123 
Pr> = J45 (262,456?, 123 + 2e2,i23e?,i2 3 — el,4 5 ei,i23 

07 = J45 



[3, 1, 1] 



— 2e|, 12361, 123 — 4e2, 4562, 12361, 123 + 4e2, 12363, 123) 



^ 2ef,i2 3 + 7e2,45el,i2 3 - 2e2, 1236^,123 + 20e 3 , me?, 123 
+10e2,i23 e i,i23 - 4e 2 , 4562, 12361, 123 - 50e 2 ,45e3,i23ei, 



12e 



2,4561,123 



15e2, 12363, 123ei, 123 



y +3ej, 45 - 10ej,i23 ~ 10e2,45el,i 23 + 15el, 123 + 18ej, 45 e 2 ,i23 
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TABLE III: Continued 
(c) N = 5 (cont.) 



A 


'-'boson 


Primitive Polynomials 


[2,2,1] 




. / ei, 3 4ef il2 - 6e?, 34 ei,i 2 + 6e 2 ,i 2 e?, 12 + 4e 2 ,34e ljl2 - 3e 134 ei,i 2 \ 
\ +13ei,34e 2 ,i2ei,i 2 + 13ei, 3462,3461, 12 + 4e 134 e 2i i 2 + 6ei, 34 e2,34 + 8e 2 ,i 2 e 2 ,34 / 
^ ei,34ei )12 + 5e? i34 e'j i il2 - 2e 2 ,i 2 ef jl2 - 6e 2 , 34 e?,i 2 + 5e?, 34 ei,i 2 - 6ei, 34 e 2 ,i 2 e?, 12 \ 

/3s = — 9ei,34e 2 ,34ef 12 + ef, 34 ei,i 2 — 4e| 12 ei,i 2 — 6el, 34 ei,i 2 — 9ei, 34 e 2 ,i2ei,i 2 — 6ei, 34 e2,34ei,i2 
\ +6e 2 ,i 2 e 2i 3 4 ei > i 2 — 6ei,34e 212 — 4ei,34e|, 34 — 6ef, 34 e 2 ,i 2 — 2ei, 34 e 2 , 34 + 6ei,3 4 e 2 ,i 2 e 2 ,34 y 
' 2e 2 ,3 4 ef il2 - ei,34e 2 ,i 2 ei jl2 + 4ei, 34 e 2 ,34e? , 12 - 2e|, 34 ei,i 2 + ei, 34 e 2 ,i 2 ei,i 2 + ei, 34 e 2 , 34 ei,i 2 \ 

A> = — 10e 2 ,i 2 e 2 , 34 ef 12 — 2ei, 3 4e|,i 2 ei,i 2 — 2ei, 3 4e|, 34 ei,i 2 + 4ef , 34 e 2 ,i 2 ei,i 2 — ei, 34 e 2 ,3 4 ei,i 2 

2 2 2 4 2 2 

y — 18ei,34e 2 ,i 2 e 2 ,34ei,i 2 — 2e 134 e 2:12 — 4e 2 ,i 2 e 2 , 34 + 2ei, 34 e 2 ,i 2 — 4e 2 ,i 2 e 2 , 3 4 — 10ei, 34 e 2 ,i 2 e 2 ,34 J 
( -3ef i34 ei,i 2 + 6e 2 ,3 4 ei,i 2 - 9ei, 34 ei,i 2 + 19ei, 34 e 2 ,34ef a2 - 9el, 34 ei il2 + 6e|,i 2 e?,i 2 + 6e|, 34 e?,i 2 \ 

2323 352 22 

+7ei i34 e 2 ,i 2 ei > i 2 + 20ei, 34 e 2 , 34 ei,i 2 — 14e 2 ,i 2 e 2 , 3 4ei,i 2 — 3e li34 ei i i 2 + 16ei,3 4 e 2 ,i 2 ei,i 2 

/Q 223232 222 

P7 = +13ei,34e 2i34 ei i i 2 + 20e 134 e 2j i 2 ei : i 2 + 7e 134 e 2 ,3 4 ei i i 2 — 42ei,34e 2 ,i 2 e 2 , 3 4ei,i 2 + 13e 134 e 2i i 2 ei,i 2 

2 2 2 4 2 2 

+ 16e li34 e 2i34 ei,i 2 — 6e 2 ,i 2 e 2 , 34 ei,i 2 + 19e lj34 e 2 ,i 2 ei i i 2 + 2e 2 ,i 2 e 2 , 34 ei,i 2 — 42e 1>34 e 2j i 2 e 2 , 34 ei,i 2 

3232 25 2 3 

\ +6ei i34 e 2i i 2 + 6ei i34 e 2i 34 + 2ei,34e 2 ,i 2 e 2 , 34 + 6ei, 34 e 2 ,i 2 — 6ei, 3 4e 2 ,i 2 e 2 , 3 4 — 14ei, 34 e 2 ,i 2 e 2 , 3 4 / 
' 3ei, 34 el jl2 + 13ef i34 ef i i 2 - 6e 2 ,i 2 ef jl2 - 22e 2 , 34 ef >12 + 30ei, 34 ei,i 2 \ 
-23ei i34 e 2 ,i 2 ei i i 2 - 84ei i 3 4 e 2 , 34 ef ji2 + 40e lj34 ef il2 + 10e| il2 ef a2 
+80e| i34 el 12 — 81e? >34 e 2i i 2 ef jl2 — 145e 1;34 e 2j34 ef , i2 + 138e 2 ,i 2 e 2 ,34ef,i 2 
+30ef i34 e?, 12 + 49ei,34e!,i 2 ef il2 + 160ei, 34 e| >34 ef , 12 

-150e? i34 e 2 ,i 2 ef a2 - 150e li34 e 2 , 34 e? >12 + 360ei,3 4 e 2 ,i 2 e 2 , 34 ef jl2 + 13ef j34 e? >12 - 10e 2>12 e ljl2 

32 222 222 22 4 2 

-22e 2 34 ei 12 + 114ei 34 e 2 i2 ei 12 + 114ei 34 e 2 34 ei, i2 - 130e 2 ,i 2 e 2 34 ei 12 - 145ei 34 e 2 ,i 2 ei i2 

4222 2 27 3 

— 81ei i34 e 2j 3 4 ei i i 2 — 72e 212 e 2j3 4ei : i 2 + 440ei, 34 e 2 ,i 2 e 2 , 3 4ei,i 2 + 3ei, 34 ei,i 2 — 37ei, 3 4e 2 ,i 2 ei,i 2 

3 3 2 3 2 2 

— 37ei,34e 2i34 ei,i 2 + 160ei i34 e 2i i 2 ei,i 2 + 49ei, 34 e 2 , 34 ei,i 2 — 270ei, 34 e 2 ,i 2 e 2 , 34 ei,i 2 
— 84ef i34 e 2 ,i 2 ei,i 2 — 23e?, 34 e 2 , 34 ei,i 2 — 270ei, 3 4el,i 2 e 2 , 3 4ei,i 2 + 360e lj34 e 2 ,i 2 e 2 , 34 ei,i 2 
-22e? i 3 4 e| il2 - 10e?, 34 e|, 34 - 44e 2 ,i 2 el, 34 + 80ef i34 e!i 2 + 10ef i34 e|, 34 - 160e| il2 el, 34 
— 72ef >34 e 2 ,i 2 e| > 34 — 22ei, 34 e 2 ,i 2 — 6ef, 34 e 2 , 34 
\ -44e 2j i 2 e 2 , 34 - 130e? i34 e| il2 e 2 ,3 4 + 138ef, 34 e 2 ,i 2 e 2 , 3 4 / 


[2,1,1,1] 




/?6 = J345 (2ef jl2 — 5e 2j i 2 ei i i 2 + e2,345ei,i2 — 2e3,34 5 ) 

pV = J345 (2e 2 ,i2ei i i 2 — e3.345e1.12 — 2e 2 . i2 + 2e 2 . 1262, 345) 

/ 2ef a2 - 7e 2 .i 2 ef. 12 — e 2 ,345ei,i 2 - 7e3,345ef,i 2 + 20e| i i 2 ei,i 2 - e|. 345 ei,i 2 \ 

P8 = J345 

\ +9e2,i2e2,345ei,i2 + 12e2, 1263,345 + 2e2, 34563, 345 / 

/ 8e 2 ,i2ei,i 2 + 2e 2 ,345ei, 12 - 863,3456?, i2 — 10e 2 ,i 2 ei, 12 + e 2 , 345 ei,i 2 — e 2 .i2e2.345e?,i 2 \ 

P9 = J345 3 2 2 2 

\ +17e2,i2e3,345ei,i2 + 4e2, 34563, 34561,12 + 20e 2 ,i 2 — 10e 2 ,i2e 2 ,3 45 — 2e 3 345 — 10e 2 ,i 2 e 2 , 3 45 / 


[1,1,1,1,1] 




J12345 
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TABLE IV: The minimal set of primitive generalized spin wavefunctions specifying a complete basis, classified by 
symmetry type A. In order to count the number of states, we count the number of primitive generalized spin wave- 
functions in each symmetry type, multiplied by the dimension of the corresponding symmetric group representation 
(see appendix [Cl for details). The total number should be equal to n N for an n component system of TV particles. 



(a) 3-Component TV = 2 



A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


[2] 


|aa> , |a/3) , |a 7 ) , |/?/3) , , | 77 ) 


6 


[1,1] 


|a/3) , |a 7 ) , |/3 7 ) 


3 



(b) 3-Component N = 3 



A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


[3] 


\aaa) , \aa/3) , |aa 7 ) , | a/3/3) , |a/3 7 ) , |a 77 ) , |/3/3/3) , |/3/3 7 ) , |/3 77 ) , | 777 ) 


10 


[2,1] 


\aa/3) , \aa-y) , | a/3/3) , |a/3 7 ) , |a 7 /3) , |a 77 ) , |/3/3 7 ) , |0 77 ) 


16 = 8 x 2 


[1,1,1] 


|a/3 7 ) 


1 


(c) 3-Component TV = 4 


A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


[4] 


|aaaa) , \aaa/3) , aaa 7 ) , aa/3/3) , |aa/3 7 ) , |aa 77 ) , a/3/3/3) , 
\aPPl) , |a/377> , l«777> , 1/5/3/3/3) , \&$M) , 1/3077) , 1/3777) , 17777) 


15 


[3,1] 


|aaa/3) , |aaa 7 ) , |aa/3/3) , |aa/3 7 ) , |aa 7 /3) , |aa 77 ) , | a/3/3/3) , |a/3/3 7 ) , 
|a0 7 0) , |a/3 77 ) , |a 77 /3) , |a 777 ) , |000 7 ) , |00 77 ) , |/3 777 ) 


45 = 15 x 3 


[2,2] 


|aa00) , |aa/3 7 ) , |aa 77 ) , |a/3/3 7 ) , |a/3 77 ) , |/3/3 77 ) 


12 = 6 x 2 


[2, 1, 1] 


|aa/3 7 ) , |a/3/3 7 ) , |a/3 77 ) 


9 = 3x3 



(d) 3-Component TV = 5 



A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


[5] 


\aaaaa) , \aaaaj3) , \aaaoi^) , aaa/3/3) , |aaa/3 7 ) , |aaa 77 ) , |aa/3/3/3) , 
|aa/3/3 7 ) , |aa/3 77 ) , |aa 777 ) , |a0000) , |a000 7 ) , |a00 77 ) , la/3777) , 
|a 7777 ) , 1/3/3/3/3/3) , |/3/3/3/3 7 ) , |/3/3/3 77 ) , |/3/3 777 ) , |/3 7777 ) , | 77777 ) 
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[4, 1] 


|aaaa/3) , \aaaa"f) , |aaa/3/3) , aaa0 7 ) , \aawyf3) , |aaa 77 ) , |aa000) , |aa/3/3 7 ) , 
|aa/3 7 /3) , |aa/3 77 ) , |aa 77 /3) , |aa 777 ) , |a/3/3/3/3) , |a/3/3/3 7 ) , |a/3/3 7 /3) , |a/3/3 7 7) , 
|a/3 7 7/3) , la/3777) , |a 7 770) , |« 7777 ) , 100007) , |/3/3/377> , l/3/3777> , 1/37777) 


96 = 24 x 4 


[3,2] 


|aaa/3/3) , |aaa07) , |aaa 77 ) , |aa/3/3/3) , |aa/3/3 7 ) , aa/377) , aa 7 00) , |aa 7 /3 7 ) , 
|aa 777 ) , |a/3/3/3 7 ) , |a/3/3 77 ) , |a/3 7 /3 7 ) , |a/3 777 ) , |/3/3/3 77 ) , |/3/3 777 ) 


75 = 15 x 5 


[3, 1, 1] 


|aaa/3 7 ) , |aa/3/3 7 ) , |aa 7 7 ) , |a000 7 ) , |a0 7 7 ) , |a 7 7/3 7 ) 


36 = 6 x 6 


[2,2,1] 


|aa00 7 ) , |aa0 77 ) , |a00 7 7) 


15 = 3 x 5 



(e) 4-Componcnt TV = 2 



A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


[2] 


|aa) , |a/3) , |a 7 ) , |a5) , |/3/3) , |/3 7 ) , |05) | 77 ) , | 7 5) , |55) 


10 


[1,1] 


|a/3) , |a 7 ) , \aS) , |/3 7 ) , \/35) | 7 5) 


6 



(f) 4-Component TV = 3 



A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


[3] 


|aaa) , |aa/3) , |aa 7 ) , \aaS) , |a/3/3) , |a/3 7 ) , \a/36) , |a 77 ) , |a 7 5) , \aS5) , 
1/3/3/3) , 1/907) , 1/3/35) , |/3 77 ) , |0 7 5) , |055) , | 777 ) , | 77 5) , | 7 55) , |555) 


20 


[2,1] 


|aa/3) , |aa 7 ) , \aa6) , | a/3/3) , |a/3 7 ) , |a05) , |a 7 /3) , |a 77 ) , |a 7 5) , |a50) , 
|a5 7 ) , \a6S) , 1/3/37) , 1/3/35) , 1/377) , l/3 7 5) , |0<5 7 ) , 1/355) , | 77 5) , | 7 55) 


40 = 20 x 2 


[1,1,1] 


|a/3 7 ) , |a/35) , |a 7 5) , |/3 7 5) 


4 
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TABLE IV: Continued 



(g) 4-Component N = 4 


A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


[4] 


|aaaa) , aaa/3) , \aawy) , \aaa8) , aa/3/3) , |aa/3 7 ) , aa/35) , |aa 77 ) , aa 7 5) , \aa58) , \aPPP) , 
loppy) , |a/3/35) , |a/3 77 ) , |a/3 7 5) , \a/366) , |q 7 7 7 ) , |a 77 5) , |a 7 55) , \aS8S) , 1/3/3/3/?) , |/3/3/3 7 ) , \P00S) , 
|£/3 77 ) , 1/9/375) , \PP66) , |/3 777 ) , |/3 77 5) , |/3 7 55) , \p5SS) , | 7777 ) , | 777 5) , | 77 55) , | 7 555) , \8585) 


35 


[3, 1] 


\aaaj3) , |aaa 7 ) , \aaa8) , aa/3/3) , |aa/3 7 ) , \aa/35) , |aa 7 /3) , |aa 77 ) , |aa 7 5) , \aa8j3) , 
\aaS-y) , \aa88) , |a/3/3/3) , |a/3/3 7 ) , \app5) , |a/3 7 /3) , |a/3 77 ) , |a/3 7 5) , |a/35/3) , |a/35 7 ) , 
\af3SS) , |a 77 /3) , |a 777 ) , a 77 5) , |a 7 5/3) , |a 7 5 7 ) , |a 7 55) , \aSS/3) , a55 7 ) , |ce5<5<5) , 
|/3/3/3 7 ) , |/3/3/35) , |/3/3 77 ) , , |/3/35 7 ) , \ppS5) , |/3 777 ) , |/3 77 5) , |/3 7 5 7 ) , |/3 7 55) , 
\pS8j) , \I35S5} , | 777 5) , | 77 55) , | 7 555) 


135 = 45 x 3 


[2,2] 


|aa/3/3) , aa/3 7 ) , |aa/35) , |aa 77 ) , |aa 7 5) , \aaS8) , |a/3/3 7 ) , |a/3/35) , a/3 77 ) , a/3 7 5) , 
|a/355) , \ajp5) , |a 77 5) , |a 7 55) , |/3/3 77 ) , |/3/3 7 <5) , |/3/355) , |/3 77 5) , \p 7 SS) , | 77 55) 


40 = 20 x 2 


[2, 1, 1] 


|aa/3 7 ) , |aa/35) , |aa 7 <5) , |a/3/3 7 ) , |a/3/3<5) , |a/3 7 5) , |a 7 /3 7 ) , |a 7 /3<5) , |a 77 5) , |a5/3 7 ) , 
|a5/35) , |a5 7 5) , |/3/3 7 5) , |/3 77 5) , |/35 7 5) 


45 = 15 x 3 


[1,1,1,1] 


|a/3 7 5) 


1 


(h) 4-Componcnt N = 5 


A 


Set of Primitive Generalized Spin Wavefunctions 


Counting 


L°J 


|aaaaa) , |aaaa/3) , |aaaa 7 ) , \aaaa8) , |aaa/3/3) , |aaa/3 7 ) , |aaa/35) , |aaa 77 ) , |aaa 7 5) , 
\aaa88) , \aa/3/3/3) , aa/3/3 7 ) |aa/3/35) , |aa/3 77 ) , |aa/3 7 5) , |aa/355) , |aa 777 ) , aa 77 5) , 
\aajS6) , |aaM5) , |a/3/3/3/3) , |a/3/3, /3 7 ) , |a/3/3/35) , |a/3/3 77 ) , |a/3/3 7 5) , \appS6) , |a/3 777 ) , 
|a, /3, 77 <5) , |a/3 7 55) , \a/3SSS) , a 7777 ) , |a 777 5) , |a 77 55) , \<yy8S5) , \aSS58) , \ppppp) , 
|/3/3/3/3 7 ) , 1/3/3/3/3,5) , |/3/3/3 77 ) , |/3/3/3 7 5) , |/3/3/355) , |/3/3 777 ) , |/3/3 77 5) , |/3/3 7 55) , \/3/36S6) , |/3 7777 ) , 
|/3 777 5) , |/3 77 55) , |/3 7 555) , \f35S58) , 77777 ) , 77 7 7 5) , | 777 55) , 77 555) , | 7 MM) , \SSS8S) 


ou 


[4, 1] 


\aaaaf3) , aaaa 7 ) , \aaaa5) , aaa/3/3) , aaa/3 7 ) , |aaa/35) , |aaa 7 /3) , |aaa 77 ) , |aQa 7 5) , \aaaS/3) , 
\aaa8^) , \aaaS8) , |aa/3/3/3) , |ao:/3/3 7 ) , \aa/3j38) , |aa/3 7 /3) , |aa/3 77 ) , |«a/3 7 (S) , |Qa/3<5/3) , |aa/3<5 7 ) , 
|aa/3(5<5) , |aa 77 /3) , |aa 777 ) , |aa 77 5) , |aa 7 <5/3) , |aa 7 <5 7 ) , \awy8S) , |aa<5(5/3) , |aa55 7 ) , \aa8S8) , 
| a/3/3/3/3) , |a/3/3/3 7 ) , | a/3/3/35) , |a/3/3 7 /3) , |a/3/3 77 ) , |a/3/3 7 5) , | a/3/35/3) , |a/3/35 7 ) , | a/3/355) , |a/3 77 /3) , 
|a,8 777 ) , |a/3 77 5) , |a/3 7 5/3) , |a/3 7 5 7 ) , |a/3 7 55) , |a/355/3) , |a/355 7 ) , |a/3555) , |a 777 /3) , |a 7777 ) , 
a 777 5) , |a 77 5/3) , |a 77 5 7 ) , |a 77 55) , a 7 55/3) , |a 7 55 7 ) , \aj8S8) , \a5S5/3) , |a555 7 ) , \aS585) , 
|/3/3/3/3 7 ) , |/3/3/3/35) , |/3/3/3 77 ) , |/3/3/3 7 5) , |/3/3/35 7 ) , |/3/3/355) , |/3/3 777 ) , |/3/3 77 5) , |/3/3 7 5 7 ) , |/3/3 7 55) , 
|/3/355 7 ) , 1/3/3555) , |/3 7777 ) , |/3 777 5) , |/3 77 5 7 ) , |/3 77 55) , |/3 7 55 7 ) , |/3 7 555) , l^^^^-y) , \I35S5S) , 
| 7777 5) , | 777 55) , | 77 555) , I'jSSSS) 


336 = 84 x 4 


[3,2] 


|aaa/3/3) , |aaa/3 7 ) , |aaa/35) , |aaa 77 ) , \aaa"/8) , \aaa8S) , |aa/3/3/3) , |aa/3/3 7 ) , |aa/3/35) , |aa/3 77 ) , 
|aa/3 7 5) , |aa/355) , |aa 7 /3/3) , |aa 7 /3 7 ) , |aa 7 /35) , aa 777 ) , |aa 77 5) , aa 7 55) , aa5/3/3) , aa5/3 7 ) , 
\aaS/3S) , \aaSjj) , |aa5 7 5) , \aaS8S) , ja/3/3^ 7 ) , |a/3/3/35) , |a/3/3 77 ) , |a/3/3 7 5) , |a/3/355) , |a/3 7 /3 7 ) , 
\af3jl35) , |a^ 777 ) , |a/3 77 5) , |a/3 7 55) , |a/35/3 7 ) , |a^5^5) , |a^5 77 ) , |a/35 7 5) , \a/35S8) , |a 77 /35) , 
|a 777 5) , |a 77 55) , |a 7 5/35) , |a 7 5 7 5) , \ ai 858) , |/3/3/3 77 ) , |/3/3/3 7 5) , |/3/3/355) , |/3/3 777 ) , |/3/3 77 5) , 
\PM58) , |/3/35 77 ) , |/3/35 7 5) , |/3/3555) , |/3 777 5) , |/3 77 55) , |/3 7 5 7 5) , |/3 7 555) , | 777 55) , | 77 555) 


300 = 60 x 5 


[3, 1, 1] 


|aaa/3 7 ) , \aaa/3S) , |aaa 7 5) , |aa/3/3 7 ) , |aa/3/35) , |aa/3 7 5) , |aa 7 /3 7 ) , |aa 7 /35) , aa 77 5) , aa5/3 7 ) , 
|aa5/35) , \aaSj5) , |a^/3/3 7 ) , |a/3^/35) , |a/3/3 7 5) , |a/3 7 ^ 7 ) , |a/3 7 /35) , |a/3 77 5) , |a^5^ 7 ) , |a/35/35) , 
|a/35 7 5) , |a 77 /3 7 ) , |a 77 /35) , |a 777 5) , |a 7 5/3 7 ) , |a 7 5/35) , |a 7 5 7 5) , |a55/3 7 ) , |a55/35) , a55 7 5) , 
|/3/3/3 7 5) , |/3/3 77 5) , |/3/35 7 5) , |/3 777 5) , |/3 7 5 7 5) , |/355 7 5) 


216 = 36 x 6 


[2, 2, 1] 


|aa/3/3 7 ) , |aa/3/35) , aa/3 77 ) , aa/3 7 5) , |aa/35 7 ) , |aa/355) , |aa 77 5) , aa 7 55) , |a/3/3 77 ) , a/3/3 7 5) , 
|a^,55 7 ) , \aPP56) , |a/3 77 5) , |a^ 7 55) , |a 7 /35 7 ) , |a 7l S55) , |a 77 55) , |^/3 77 5) , |/3/3 7 55) , |/3 77 55) 


100 = 20 x 5 


[2,1,1,1] 


|aa/3 7 5) , |a/3/3 7 5) , |a 7 /3 7 5) , a5/3 7 5) 


16 = 4 x 4 
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TABLE V: Linear combinations of primitive 3-componcnt generalized spin wavefunctions that lead to spin-1 eigen- 
functions. In the spin language the notation should be interpreted as: a, /3 and 7 being s z = +1, and —1 respectively. 
Note that only highest possible S z states of a given S series are listed: additional states with the same S eigenvalue 
can be found by applying a spin lowering operator to any of the states listed here. Additionally we have listed the 
total number of states for each symmetry type (including the S z degeneracy): the total should equal the number of 
3-component generalized spin wavefunctions for that symmetry type given in table IIVI 



(a) N = 3 





A 




[3] 


[2,1] 


[1,1,1] 


3 


\aaa) 






2 




\aaP) 




1 


| a/3/3) - 2 \aa-y) 


2 \aPP) - \aaj) 











\apj) 


Allowed S 


1,3 


1,2 





Counting 


10 = 7 + 3 


8 = 5 + 3 


1 



(b) N = 4 





A 


s z 


[4] 


[3,1] 


[2,2] 


[2, 1, 1] 


4 


\aaaa) 








3 




\aaaP) 






2 


\aaPP) — 2 |aaa7) 


| aaPP) — | \aaa^) 


\aapp) 




1 




4 \aPPP) - \aaPf) - 3 \acryP) 




\aaP~y) 





4 |aa 7 7> - 4 \app-y) + 1/3/3/3/3) 




4 \aPP-i) — aa77) 




Allowed S 


0,2,4 


1,2,3 


0,2 


1 


Counting 


15 = 9 + 5 + 1 


15 = 7 + 5 + 3 


6 = 5 + 1 


3 



(c) N = 5 





A 


s z 


[5] 


[4,1] 


[3,2] 


[3, 1, 1] 


[2, 2, 1] 


5 
4 
3 

2 

1 



\aaaaa) 
\aaaPP) — 2 \aaaaj) 

4 1 00077) — 4 \ aaPP'y) 
+ \aPPPP) 


|oaoo/3) 
\aaaPP) — | |aaao7) 
\aaPPP) - i |oaa/37) 
-if \aaa-yP) 
\aPPPP) - \aap 1 p) 
— 5 |aa/3/37) + | |qqq77) 


|aaa/3/3) 

\aaPPP) - | |aaa/37) 

|aa7/3/3) + | |aaa77) 
— |aa/3/37) 


\aaap-y) 
a/3/3/37) — 007/37) 


\aappj) 


Allowed S 


1,3,5 


1,2,3,4 


1,2,3 


0,2 


1 


Counting 


21 = 11 + 7 + 3 


24 = 9 + 7 + 5 + 3 


15 = 7 + 5 + 3 


6 = 5 + 1 


3 
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TABLE VI: Decomposition of SU(4) appropriate for the cases such as graphenc which has SU(2) x SU(2) internal 
states (spin x valley). In the spin language the notation should be interpreted as a is s z = +1/2 in first valley and /3 
is s z = —1/2 in first valley; 7 is s z = +1/2 in the second valley and 5 is s z — —1/2 in the second valley . Note that 
only the highest possible S z states of a given S are listed: additional states with the same S eigenvalue can be found 
by applying a spin lowering operator to any of the states listed here. We have also listed the total number of states 
for each symmetry type (including the S z degeneracy): the total should equal the number of 4-component generalized 
spin wavefunctions for that symmetry type given in table IIVI 





(a) JV 


= 2 




A 


s 


[2] 


[1,1] 


1 


I77) , \aj) , \aa) 


\ori) 





\Pi) - \aS) 


| 7 *> , \M) ~ M) , \afi) 


Total No. 


10 = 3 x 3 + 1 


6=1x3+3x1 



(b) N = 3 





A 


s 


[3] 


[2,1] 


[1, 1, 1] 


3/2 


1777} , l"77> . \aary) , \aaa) 


|a77) , |aa7) 




1/2 


\aa5) + 2 |a/37) , \ aa P) 


|77<5) , ^77) - |<*7<S) , \a6-y) - \a-yS) 
2 \aaS) + \oeyf3) , \af3j) — \aa5) , aa/3) 


\a-yS) , |a/37) 


Counting 


20 = 4x4 + 2x2 


20 = 2x4 + 6x2 


4 = 2x2 



(c) N = 4 





A 


s 


[4] 


[3,1] 


[2,2] 


[2, 1, 1] 


[1,1,1,1] 


2 


]777 7 ) , 10777) , 
|aa77) , \aaaj) , 
|aaaa) 


l Q 777) , 1^077) , |aaa7) 


|aa77) 






1 


1^777) - \a-yyS) , 
la/377) — , 
laa/37) — |aaa<5) 


|777<5) , |aaa/3) , 
l^777> - |as77<5) , 
|a7^7) — \ofyy8) , 
\awyS) + 077,8) , 
la/577) — \ct-Ct-y5) , 
| aa<57) — \aotry8) , 
| \aaa8) + \aaj/3) , 
laa/37) — \aaaS) 


l«77^) , 

la/377) — |aa7<5) , 
laa/37) 


|o77<5> , 
1 07^7) , 
|aa7<5) , 
aa/37) 







I/3/377) - 2 \afijS) 


\eeyS5) — \aSSj) 

- |/3 77 5) + IfaSj) , 
\aa5S) — | a/?75) 

- |a/3<5 7 ) + I/3/377) , 
-1 \aaf38) - | \aa5f3) 
+1 |a/3/3 7 ) + |a/3 7 /3) 


I77M) , |aa/3/3) , 
1/877*) - I \arySS) , 
\aaSS) - 2 |a/3 7 (5) + l/J/377) , 
1 \aPj6) + |a 7 /3<5) , 
|a/3/37) - \ \aa0S) 


|/?775) — laS'yS) , 

\aSf3-y) - |a/3 7 5) - \ajf3S) , 

|a/3/37) - \ \<*aPS) 


|a/3 7 <5) 


Count. 


35 = 5x5 + 3x3 + 1 


45 = 3x5 + 9x3 + 3 


20 = 1x5 + 3x3 + 6 


15 = 4 x 3 + 3 


1 



(d) iV = 5. In this table we simply state the 
number of cases. 





A 


s 


[5] 


[4, 1] 


[3,2] 


[3, 1, 1] 


[2, 2, 1] 


[2,1,1,1] 


5/2 


6 


4 


2 











3/2 


4 


12 


6 


6 


2 





1/2 


2 


6 


12 


6 


6 


2 


Counting 


56 


84 


60 


36 


20 


4 
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IV. DISCUSSION 

The main result of this work is the enumeration of the 
space of TV-particle wavefunctions for bosons or fermions 
with internal degrees of freedom such as spin (i.e., multi- 
ple components). These wavefunctions can be used as a 
basis for writing pseudopotential coefficients to parame- 
terize physical problems, as, for example in Refs. and 
l4rl More importantly, however, determining such a com- 
plete basis strongly suggests much of the physics of wave- 
functions that can be generated with model interactions. 

To elaborate on this last point, let us return to the 
spinlcss (single component) case for a moment, and for 
simplicity, let us consider bosons. In the case of TV = 3 
body interaction, there are 3-body wavefunctions with 
L = 0, 2, 3, 4 . . . (this can be read off of the table from 
Ref. [l3|, which is identical to the TV = 3 fully symmet- 
ric A = [3] case of table U in the current paper). One 
can consider a family of model Hamiltonians that succes- 
sively gives each of these terms some positive energy^ 7 -. 
For example, if L = is given positive energy (i.e., there 
is a nonzero pseudopotential coefficient Vl=o,n=3) and 
no other potential energy term is in the Hamiltonian, 
then the ground state is the Moore- Read Pfaffian&, which 
has the property that when three particles approach each 
other, the wavefunction vanishes as L = 2 powers (since 
L = 1 is forbidden) . One can then consider adding a pos- 
itive energy for L = 2 as well (a positive pseudopotential 
Vl = 2.n=3), which results in the Gaffnianii as its ground 
state wavefunction, where the wavefunction vanishes as 
3 powers when 3 particles approach each other. One can 
continue by adding positive energy for L = 3, resulting 
in the Haffnianis which then vanishes as 4 powers when 
three particles approach each other. For higher L the 
situation becomes somewhat more complicate d 15 ' 17 for 
several reasons. First, there can be several wavefunctions 
with the same value of L (e.g., two such TV = 3 wavefunc- 
tions at L = 6), so that specifying L alone (equivalent to 
specifying a thin-torus limi t 55 ' 56 or a pattern of zeros 5 ^) 
does not fully specify the wavefunction. 

Secondly, it is possible that a Hamiltonian which for- 
bids all L < Lq for clusters of TV particles can have a 
ground state wavefunction where L > Lq for clusters of 
TV particles rather than L = Lq (see for example the dis- 
cussion in Ref. [j~7h . Nonetheless, the general idea that 
one can dictate the vanishing properties of a wavefunc- 
tion by appropriately choosing nonzero pseudopotentials 
remains a powerful approach both to understanding the 
properties of quantum Hall wavefunctions, and also to 
generating new and interesting trial states. 

To generalize this idea to the multi-component case, 
one can first look at table fl] If one is interested in 
two component (spin-1/2) quantum Hall states, for say 
bosons, we can look at table Hlb and quickly see what 
kind of pseudopotentials are possible — which then also 
suggests what kind of wavefunctions might occur. In the 
simplest case we might consider two-body interactions. 
Trivially in this case (looking in the TV = 2 rows of the 



table) we see that even L occurs in the triplet channel 
(S = 1) whereas odd L occurs in the singlet channel 
(5 = 0). Choosing to give positive energy (positive pseu- 
dopotential coefficient) to all L < m (triplet) even and all 
L < n (singlet) odd generates the Halperin mmn ground 
stated. 

Let us now move on to three-body interactions and 
for simplicity let us still consider spin-1/2 bosons. In the 
spin polarized (5 = 3/2) channel we see exactly the same 
structure we did for spinless bosons — wavefunctions oc- 
curring at L = 0, 2, 3, 4 . . . (TV = 3, 5 = 3/2 line of table 
HJb). In the 5 = 1/2 channel on the other hand, we see 
wavefunctions at L = 1,2,3.... A simple example of a 
three-body Hamiltonian is one that forbids L = in the 
polarized (5 = 3/2) channel (a positive V^~ Q 3 ^ =3 only); 
such a Hamiltonian generates the k = 2 NASS stat e - 38 ' 39 
as its ground state. As one might suspect, this wavefunc- 
tion vanishes as L = 2 power when three particles come 
together in the 5 = 3/2 channel and vanishes as L = 1 
power in the 5=1/2 channel. These are the lowest pow- 
ers on the table that are not explicitly forbidden by the 
Hamiltonian. 

By similarly examining the tables presented in the cur- 
rent work we can easily propose new wavefunctions that 
generalize those already discussed, and we may even be 
able to guess at the model Hamiltonians that generate 
these wavefunctions as their ground states. For exam- 
ple, generalizing the k = 2 NASS state we might pro- 
pose to forbid the lowest powers that are not forbidden 
by the NASS hamiltonian — i.e., we choose to forbid 
TV = 3,5 = 3/2, L = 2 and the TV = 3, 5 = 1/2, L = 1 
in addition to the TV = 3, 5 = 3/2, L = of the NASS 
wavefunction. We might guess that the resulting wave- 
function should vanish L = 3 powers when three particles 
come together in the 5 = 3/2 channel and as L = 2 pow- 
ers when three come together in the 5=1/2 channel. 
Such a wavefunction, a spin-singlet generalisation of the 
Gaff man, has been discussed recently in Refs. |44| andl58l. 

The examples discussed here are only a fraction of 
the model Hamiltonians that have been discussed in 
the literature^^ i 32 1 35 1 38-40 1 4ai43. Since the language of 
multicomponent pseudopotentials is arbitrarily general, 
all of these model Hamiltonians can be rephrased into 
this language. In many cases it is quite obvious from 
looking at our tables what vanishing behaviors for TV- 
particle cluster are being forbidden and what the result- 
ing properties of the corresponding ground state wave- 
functions should therefore be (albeit rigorous proofs of 
these correspondences may be more tricky). 

Once one starts considering multicomponent wavefunc- 
tions, there are clearly quite a few model Hamiltonians 
and correspondin g m odel wavefunctions that can be con- 
sidered (see Refs. l44l and[45| for recent progress in these 
directions). For example, in Ref. spin-singlet versions 
of Jack wavefunctions have been constructed. However, 
a much larger variety of states may be constructed by 
relaxing some of the restrictions of that work. One pos- 
sible direction is the generation of wavefunctions which 
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mix N and N' body multiconiponcnt interactions — for 
example, the double PfafRan wavefunction discussed in 
Ref. |32 is the ground state of the sum of a three body 
interaction and a two body interaction. Because of the 
added richness of the multicomponent case over the sin- 
gle component case there are certainly a far wider variety 
of possible model Hamiltonians and model wavefunctions 
to be explored. Potentially we may even find some new 
physically realizable multicomponent wavefunctions with 
interesting braiding properties which can be exploited for 
quantum information processing^. 

Even in the simple spinless (single component) case 
there remain quite a few open questions about model 
wavefunctions, such as Jacks, and model Hamiltonians; 
and these questions generalize to similar questions about 
the multicomponent cases. For example, it remains un- 
certain in many cases whether simple (physical) Hamilto- 
nians exist that generate these special wavefunctions as 
their ground states. Further, outside of the simplest Jack 
wavefunctions (the Read-Rezayi series and the Laughlin 
wavefunctions) , all of the spin-polarized Jacks correspond 
to nonunitary conformal field theorie o 59 i 60 which excludes 
them from describing gapped states of matter—. In at 
least some case s 11 ' 12 it is suspected that these special 
wavefunctions may correspond to critical points between 
gapped phases — however, a general understanding of 
the criticality of these special wavefunctions is lacking. 
There have been, however, some attempts to understand 
some of these properties for the spinless case by under- 
standing the relationship to the multicomponent cases 
first^ 2 -. This provides yet another motivation for explor- 
ing multicomponent Hamiltonians and wavefunctions. 
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APPENDICES 

In this set of appendices we shall set out to provide a 
detailed exposition of the linear vector space inhabited 
by (multi-component) Haldane pseudopotentials, that is, 
the vector space defined in Eq. Q21 

|i,q>. 

In doing so, the key realization is that vector space of 
multi-component pseudopotentials is identical to the vec- 
tor space of co-ordinate wavefunctions describing a multi- 
component quantum Hall system. We can thus map the 



problem of discovering a suitable basis in which to de- 
scribe pseudopotentials in the LLL, to a problem of con- 
structing multi-particle co-ordinate wavefunctions in the 
LLL with an internal degree of freedom. The properties 
required of these wavefunctions are that they should be 
homogeneous of degree L, translationally invariant, and 
overall (anti)symmetric for (fermions)bosons. These con- 
ditions are equivalent to the statement that each wave- 
function would be a rotationally invariant state on an 
appropriately sized sphere. 

In these appendices we shall provide a description of 
how to construct such wavefunctions ip = \L,q). Very 
generally, systems with n internal degrees of freedom are 
characterized by irreducible representations of the special 
unitary group SU(n). The fundamental objective of this 
work is thus to decompose the space of wavefunctions 
ijj in terms of a basis of wavefunctions corresponding to 
irreducible representations of SU(n). 

The layout of the appendices will be as follows: in 
appendix [X] we shall allude to some important mathe- 
matical preliminaries, which are to be employed in the 
remaining appendices; in appendix [B] we shall describe 
the general procedure for the construction of the basis 
of wavefunctions from a combination of spatial functions 
and generalized spin wavefunctions ; in appendix [C] wc 
shall concentrate on the details of the generalized spin 
wavefunctions and we shall derive the results presented in 
tables lIVl IVl and IVT1 in appendix IDl we shall concentrate 
on the details of the spatial parts of the wavefunctions, 
in particular we shall apply the procedure of appendix 
iBlin the context of quantum Hall wavefunctions, and we 
shall derive the co-ordinate wavefunctions listed in ta- 
bles [IT] and IIIII and the dimensions of the vector space 
of polynomials listed in table fl] finally, appendix lEl deals 
with some more advanced mathematical underpinnings 
of these ideas. 



Appendix A: Mathematical Prelimaries 

Our goal in this appendix will be to deal with the 
underlying mathematical concepts and theorems which 
are to be employed in later appendices. The underly- 
ing mathematics is that of integer partitions, of Young 
Tableaux and of the representation theory of the sym- 
metric group and special unitary group SU(n). 



1. Integer Partitions 

The concept of an integer partition is frequently used in 
this work. An integer partition A is defined as follows^: 
for a positive integer N, an integer partition is a way 
of writing N as a list of k positive integer summands 

A = [jVi,JV 2 ,...iVfc] for which Yn=i N i = N - The list 
of summands is typically presented in weakly descending 
order (that is N x > N 2 > N 3 . . . > N k ). For example, 
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the set of integer partitions for the integer 5 are [5], [4, 1], 
[3,2], [3,1,1], [2,2,1], [2,1,1,1], [1,1,1,1,1]. 

When writing down a set of integer partitions such as 
that given in the above example it is clearly important 
to specify an ordering system. In this paper we shall 
make use of lexicographic ordering^: when comparing 
two integer partitions A and fj, we say that A occurs before 
/2, denoted A < fi, if A, = /Xj for i = 1, for some 
integer j between 1 and k — 1, and then > Mj+i- In 
the above example we have used lexicographic ordering, 
for instance [3, 2] occurs before [3, 1,1]. We shall discuss 
ordering methods in more detail in appendix [El 

2. Young Tableaux 

The concept of Young tableaux is an important tool for 
categorizing group representations, which we shall need 
to do shortly. Before discussing specific group representa- 
tions, let us first explain how to construct Young tableaux 
in general. 

To draw a Young tableau, one first draws a Young 
frame or diagram. A Young frame is constructed as fol- 
lows: given an integer partition A, we build a frame by 
placing a series of empty boxes in rows and columns. 
There is one row for each of the k summands in the in- 
teger partition: the ith row contains Xi boxes. We shall 
use the convention that a new row is started beneath the 
previous row and that the rows are in the order speci- 
fied by the integer partition, hence the length of a row 
will always be less than or equal to the length of the row 
above. In addition, the rows are to be left-justified. For 
example, the Young frame corresponding to the integer 
partition A = [3, 2, 1, 1] would be: 



A Young tableaux can be constructed from a Young 
frame by placing a series of ./V integers in the boxes. 
The set of integers can be chosen in several ways — of 
particular significance to our work are the arrangements 
called standard tableaux and the arrangements called 
semi-standard tableaux. 

In a standard tableaux, the integers 1 to N are placed 
in the boxes of a Young frame of size N, with each num- 
ber occurring precisely once. Further, we place the num- 
bers such that the integer placed in any box is strictly less 
than both the integer placed in the box immediately to 
the right and the integer placed in the box immediately 
below. For example, here is a possible standard tableau 
constructed from the Young frame given above: 



1 


2 


3| 


4 


5 




6 






7 







There can of course be multiple ways to satisfy these 
conditions, in other words, if the shape of a Young 
tableaux is specified by an integer partition A then there 
is a corresponding set of size f x of admissible standard 
tableaux. We introduce an index r, which runs from 1 to 
/ , to distinguish between the possible standard tableaux 
T r in the set. 

In order to specify one particular standard tableau in 
the set of T r we must first order that set: to do so we 
shall again employ a kind of lexicographic ordering sys- 
tem. For a given standard tableau we generate a list of 
numbers starting with the number in the top left box (for 
the above tableau that would be 1), then working along 
each row from left to right taking the numbers from these 
boxes and then adjoining the numbers from each subse- 
quent row in the same fashion. Doing so will generate a 
list of N numbers Ni to Nn (in the above tableau that 
would be [1, 2, 3, 4, 5, 6, 7] ) and we then simply apply the 
procedure of lexicographic ordering to these lists. 

Let us demonstrate the definitions given in the last two 
paragraphs via the following example: the integer parti- 
tion A = [3, 1] has /I 3 ' 1 ' = 3, and the lexicographically 
ordered set of standard tableaux are 



1 


3 | 4 | 


1 


2 | 4 | 


1 


2 | 3 | 
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A semi-standard tableaux is constructed along similar 
lines: in this case however one is allowed to place any 
of the numbers 1 to n in the boxes of a Young frame 
of size N, where N does not have to equal to n. It is 
not necessary to include every number from 1 to n and 
the numbers can also be repeated. Further, we place the 
numbers such that they are non-decreasing from left to 
right and strictly increasing from top to bottom. For ex- 
ample, here is a possible semi-standard tableaux of type 
n = 4 and N = 8: 



1 


1 


4 4 


2 


2 




3 


3 



The complete set of semi-standard tableaux for a given 
A is enumerated by listing every possible arrangement of 
subsets of the numbers 1 to n in the N boxes and satisfy- 
ing the conditions described above. As with the standard 
tableaux, there are multiple ways to satisfy these condi- 
tions and we therefore have a set of h x admissible semi- 
standard tableaux. Semi-standard tableaux can also be 
ordered lexicographically, using the same procedure as 
for the standard tableaux. For example, with A = [2, 1] 
and with n = 3 we have hP 1 ' 1 ^ = 8 and the set of lexico- 
graphically ordered semi-standard tableaux are: 



1 


l| 


1 


2\ 


1 


3| 


1 


1 


1 


2 


1 


3| 


2 


2| 


2 


3| 


2 




2 




2 




3 




3 




3 




3 




3 





Finally, we shall introduce the term conjugate to refer 
to any tableaux related by a reflection along the diagonal, 
for example, the following tableau of shape A = [3, 1]: 



1 


3 | 4 | 
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is conjugate to the following tableau of shape A = [2, 1, 1] 



1 


2 


3 




4 





3. Irreducible Representations of the Symmetric 
Group 

The symmetric group Sjv is the group of permutations 
on N objects, {tti, ...,7Tjvi} • A permutation 7T; of labels 
(123. ..N) to (iii2i3— ijv) is denoted by (iii2i3---ijv)- A 
transposition (pair interchange) of label i and label j with 
all other labels unchanged is denoted by (i;j). 

The classification of multi-component wavefunctions 
depends essentially upon the properties of Sat, and in 
particular on the understanding of irreducible represen- 
tations of Sat. Throughout this appendix we shall as- 
sume that the reader is familiar with the fundamental 
concepts of representation theory; a more detailed dis- 
cussion can be found in Ref. l48l The correspondence 
between wavefunctions and representations of the sym- 
metric group manifests itself in the concept of a symme- 
try type, as described in section IIII Bl The interpreta- 
tion in terms of symmetry types is said to be the physi- 
cists' interpretation of the mathematics of the irreducible 
representations^. 

In order to construct multi-component wavefunctions, 
we will find that we need to construct an appropriate 
symmetric group algebra. A well known example of a 
symmetric group algebra is the algebra of Young op- 
erators (for their definition see appendix [E)l . There is, 
very generally, a direct correspondence between a repre- 
sentation of a group and a representation of the group 
algebra 4 ^. For our purposes it will be convenient to con- 
struct the symmetric group algebra corresponding to the 
orthogonal representation of the symmetric group (the 
reasoning behind this choice is explained in appendix 
|E| . We shall define the orthogonal representation of the 
symmetric group shortly. In order to construct such 
an algebra we will then make use of the matric basis 
of the symmetric group algebra, which is a very gen- 
eral method to construct a basis directly from a group 
representatio n 64 ! 65 . 



representation matrices are themselves orthogonal matri- 
ces, meaning that they satisfy U (tt," ) = U (7T,) sr . 

In fact, orthogonal representation matrices can be con- 
structed directly from standard tableaux using a measure 
called the axial distance^ 5 -: suppose that the number p 
appears in the i p th row and the j p th column of a given 
standard tableau and that the number q appears in the 
iqth row and the jqth column. The axial distance from 
p to q in tableaux T r is defined as 

^P,q = (%> — ip) — (*'q — Jq) ■ 

Using this definition one can construct orthogonal rep- 
resentation matrices U (n) for 7r £ Sjv as follows: first 
we define the elements of the representation matrix 
U ((k - 1; k)) of transpositions (k - 1; k) for k = 1 . . . N: 



Z7((k-l;k)) w = -l/d 



k-l,k 



tf((k-l;k)) 4 



1 



if T i9 6(k-l;k)T< 
if T 3 = (k-ljk)!}. 



Then we note that every permutation can be expressed 
as a product of transpositions of the form (k — 1; k); the 
proof is given in Ref. [65| for example. Using this result 
we can generate the representation matrices correspond- 
ing to every other element of the symmetric group by 
simple matrix multiplication. This is Young's orthogo- 
nal representation. Some of these orthogonal matrices 
are tabulated in Ref. |H. 

Representations of the symmetric group correspond- 
ing to conjugate Young tableaux shapes are known as 
contragradient representations (see Ref. We can 

explicitly construct contragradient orthogonal represen- 
tation matrices V (ni) X , which are related to the original 
orthogonal representation matrices U (iri) X by 

V{n)l s = sgn(7Ti)C/ (n)r a ■ 

Here we have used the sign of a permutation, sgn(7Ti), 
which is defined as follows: writing a permutation as a 
product of transpositions, if there are m transpositions in 
the product then the sign of the permutation is sgn (wi) = 
(-l) m 



b. The Matric Basis of the Symmetric Group Algebra 



a. Young 's Orthogonal Representation of the Symmetric 
Group 

A fundamental principle of the representation theory 
of the symmetric group is that the irreducible represen- 
tations of the symmetric group are in one-to-one corre- 
spondence with the standard tableau?* 4 ^. Given the set 
of Young tableaux of shape A we can always construct 
some kind of irreducible representation matrices U (7i"i) A 
of dimension / . In order for the representation to be an 
orthogonal representation, all that is required is that the 



For any given set of symmetric group representation 
matrices U (7i"i) A of dimension / A we can in general form 
the following elements of the corresponding symmet- 
ric group algebra, called matric units 6 ^ or semi-normal 
units^ 5 -: 

S "s = ^ E U K rl ) I *i r,s= 1, f\ (Al) 
i=l 

These elements are linearly independent (the proof of 
which is given in Ref. [65h and can be chosen as a basis 
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for the symmetric group algebra (see Ref. They 
satisfy the orthogonality condition (with no summation 
implied) : 

&rs&uv = ^X^su^rv (A2) 

The diagonal elements, e rr , are idempotent and mutually 
orthogonal. Some explicit forms of these operators are 
tabulated in tabic IIXI in appendix [E] 

Using the matric units, we can also define a basis for 
the group algebra of a cotragradient representation: 

/ A N ' A 
^» = j^_J2 V ^rs■ K ^■ 
i=l 

and these elements satisfy the same condition as above, 
namely: 

-p. — X. X pA 

4. Irreducible Representations of SU(n) 

The special unitary group, denoted SU(n), is the group 
of n x n unitary matrices with determinant 1. SU(ri) also 
describes the structure of multi-component systems. Ir- 
reducible representations of SU(?i) are in one-to-one cor- 
respondence with certain subsets of the irreducible rep- 
resentations of the symmetric group — specifically those 
representations associated with the set of standard Young 
tableau restricted to having at most n rows (see Ref. |48| 
for more details). 



5. The Lie algebra of SU(n) 

In this work we describe wavefunctions with internal 
degrees of freedom such a spin. The reader will likely 
be familiar with the fact that quantum states labeled 
by spins are classified by the irreducible representations 
of the Lie algebra of SU(2), that is the algebra of the 
S 2 and S z operators. More generally, a generalized spin 
wavefunction will be classified by representations of the 
Lie algebra of SU(n). Irreducible representations of the 
Lie algebra of SU(n) are called multiplets and each mul- 
tiplet contains a vector space of associated functions or 
states^. For example, for SU(2) the multiplets corre- 
spond to the values of the quantum numbers S and the 
states within each multiplct are labeled by the possible 
values of S z ; for SU(3) the multiplets correspond to the 
classification of baryonic particles in the quark model^. 
The general form of the SU(n) multiplets is described in 
detail in Ref. [66|. Since the multiplets correspond to the 
irreducible representations of SU(n), they are labeled in 
accordance with the set of Young tableaux of shape A. 
In other words a multiplet labeled by A corresponds to a 
particular symmetry type A. 



In order to label a particular state in a given multiplct, 
one must specify the component content of that state: for 
example, in the two-component case, the number of spin 
up and spin down particles (or alternatively the S z eigen- 
value); in the three-component case, the quark content 
of the state, and so on. We shall introduce the notation 
A z , which will specify the component content of a given 
state, and thus distinguish between different states in the 
same multiplct (the nomenclature is inspired by S z ). In 
addition, there can be multiple states which occur for the 
same value of A z — the number of such states is called 
the multiplicity M(A, A z ). 

There is a one-to-one correspondence between the 
states of a SU(n) multiplet and the set of semi-standard 
tableau:?*^. For a given A we construct the set of semi- 
standard tableaux containing up to n different numbers 
— these can be interpreted as the set of components A 2 
for that state. Recall our earlier example of the set of 
SU(3) semi-standard tableaux: 



f 


l 
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corresponds to X z = {2,1,0}, and since there is 
only one occurrence of this type if tableaux we have 
M([2,l],{2, 1,0}) = 1. Similarly we have 



1 


2 


1 


3 


3 
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which both correspond to X z = {1,1,1}. In this case 
we clearly have M([2, 1], {1, 1, 1}) = 2 for SU(3). The 
total number of states in each multiplet, including the 
multiplicity, is the SU(n) dimension of the representation 
it is calculated by enumerating the complete set of 
semi-standard Young tableaux for a given A, for example, 
with A = [2,1] and n = 3, the dimension is 8. 

A more complete description of the Lie algebras of 
SU(n) and their representations can be found in Ref. 16(1 . 

6. Tensor Products of SU(n) Multiplets 

The problem of constructing a vector space of all pos- 
sible N-particle generalized spin wavefunctions is equiva- 
lent to the mathematical problem of forming tensor prod- 
ucts of N fundamental multiplets. A fundamental multi- 
plet is the irreducible representation corresponding to a 
Young tableaux containing a single box — so for SU(ra) 
the dimension of the fundamental multiplet is equal to 
n. Such tensor products can be decomposed into a direct 
sum of multiplets^.. An important feature is that in this 
direct sum it is possible to have multiple occurrences of 
the same multiplet (we shall give an example shortly). 
In our general construction we shall therefore introduce 
an index r which will distinguish between such repeated 
multiplets. This decomposition enables us to determine 
the possible multiplets and multiplet states which fully 
span our vector space. 
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For illustrative purposes we shall go through the three- 
component case in detail here. A three-component sys- 
tem can be decomposed into SU(3) multiplets. We can 
determine the allowed set generalized spin wavefunctions 
by forming tensor products of SU(3) representations of a 
single particle state. 

We start with a single particle state with three com- 
ponents, which corresponds to the fundamental SU(3) 
multiplct. To deduce the vector space of two-particle 
generalized spin wavefunctions we form a tensor product 
of two fundamental multiplets (here we label the multi- 
plets by the shapes of the corresponding young frames): 



e 



This product is evaluated according to the Littlewood- 
Richardson rule, see for example Ref. HH to give a direct 
product of SU(3) multiplets. The [2] representation is 6 
dimensional, whilst the [1,1] representation is 3 dimen- 
sional (c.f the N = 2 case in table HV]). 

If we want to include a third particle then we construct: 



□®n®n=i 



We see that there are two copies of the [2,1] multiplct, 
which has dimension 8. The two copies are distinguished 
by the index r and the states within the two multiplets 
are mutually orthogonal to each other (so in this example 
we can have r = 1 or r = 2). The [3] multiplet has 
dimension 10, and the [1,1,1] has dimension 1: so we see 
that the total dimension is 10 + 2 x 8 + 1 = 27 = 3 3 
(c.f. the N = 3 case in table ITV|) . The total dimensions 
calculated in this way are consistent with the counting of 
states listed in table IIVI 

We can also apply the Littlewood-Richardson rule to 
the construction of standard tableaux, thus we observe: 



[T]®[2]<8 r jrp 
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1 


3 




i 


2 




1 1 2 | 3 




_2_ 




_3_ 





We see that the two copies of the [2,1] multiplet can 
be distinguished by their association to standard Young 
tableaux. This result comes about because of the un- 
derlying connection between tensor representations of 
groups and irreducible representations of the symmetric 
group: specifically the symmetry property of the indices 
of the tensor itself forms a representation of the symmet- 
ric group. 

More generally therefore, the number of repeated oc- 
currences of a particular Young tableaux shape A simply 
corresponds to the dimension f x of the symmetric group 
representation which would also be denoted by a Young 
tableaux of that shape (a proof of this result is given, for 
example, in Ref. |48|, Ch. 10). For example, the symmet- 
ric group representation dimension of the [2,1] tableaux is 



2 (recall that the symmetric group representation dimen- 
sion is given by the number of standard tableaux for a 
given tableaux shape). Thus one can associate a specific 
standard Young tableaux shape with each occurrence of 
an SU(re) multiplct in our decomposition. 

To summarize, the tensor product of N fundamental 
multiplets will contain every admissible irreducible rep- 
resentation of the Lie algebra of SU(n), that is, the irre- 
ducible representations corresponding to every possible 
standard tableaux of shape A with size N and with no 
more than n rows^, and in such a decomposition the 
multiplct described by A will occur / A times, with the oc- 
currences being distinguished by an association to stan- 
dard Young tableaux. 



Appendix B: General Procedure for Construction of 
Multi-Component Wavefunctions 

In this appendix we shall describe the solution to the 
very general problem of constructing wavefunctions for 
multi-component systems by decomposing into a basis 
of spatial wavefunctions and generalized spin wavefunc- 
tions. In the following appendix, appendix [Cj we shall 
describe the forms of the generalized spin wavefunctions 
in more detail. Finally, in appendix [D] we shall describe 
how to apply the general arguments for construction of 
spatial wavefunctions presented in the current appendix 
to the special case of construction of translationally in- 
variant LLL quantum Hall wavefunctions with fixed rel- 
ative angular momentum. 

The vector space of wavefunctions describing a system 
of identical particles with internal degrees of freedom can 
always be written as a tensor product of the spaces of 
spatial and generalized spin wavefunctions: 



spatial/ 



\x 



gsw/ 



Thus the dimension of the Hilbert space (with no symme- 
try restriction) is equal to the product of the dimensions 
of the space of spatial and generalized spin wavefunc- 
tions. Additionally such a wavefunction must conform to 
the symmetry condition required for a system of identical 
fcrmions (bosons) — namely the wavefunction must be 
antisymmetric (symmetric) under any exchanges of the 
particle labels. 

If we are trying to determine the dimension of the space 
of wavefunctions with N particles of some degree L which 
obeys the proper symmetry condition we can in principle 
start with the most general spatial function 3> sp atiai of 
this degree, and the most general generalized spin wave- 
function X gsw (where a "most general wavefunction" con- 
tains one arbitrary coefficient multipying each of the basis 
functions in the space) . The most general $ S patiai is given 
by a linear combination of all possible arrangements of 
N particles into ./V orbitala single-particle orbitals <f>\ ...4>n 
(including multiple occupancy) with a fixed combined or- 
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der: 

Norbitals 

^spatial = ]P a lu ....i N (<j)i 1 (r 1 )...(/) iN (r N )), (Bl) 

ii,...ijv = l 

where fl^ are arbitrary coefficients. Similarly the 
most general X gsw is given by 



X, 



where arbitrary coefficients and where §j are 

the complete set of generalized spin wavefunctions of 
the form given in Eq. 1231 of which there are n N for 
an iV-particle, n-componcnt system. The most general 
combined wavefunction is given by the product of the 
generalized spin wavefunction and spatial wavefunction: 

Agsw ^spatial ■ 

Next we impose that the wavefunction ip is overall anti- 
symmetric (for fermions) by applying an antisymmctrizer 
A to the most general construction 



Ip = -4.$ S patial^g 



(B2) 



where 



Nl 



A = 



^sgn(7r ? ;)7ri. 



For bosons we simply replace the antisymmetrizer by a 
symmctrizer S: 



S 



nsn. 



i=l 



After applying one of these symmetry operators to the 
most general construction we will be left with a wavefunc- 
tion containing a certian number of linearly independent 
coefficients: the total number of those coefficients is the 
size of the Hilbert space. One can always deduce by brute 
force the number of linearly independent coefficients in 
ip by determining the linearly independent combinations 
of the a,i lt ...,i N and Xj coefficients remaining after the 
symmetry operator is applied. In general, however, this 
method is extremely inefficient and becomes essentially 
algebraically impossible for more than a very few parti- 
cles and very low degrees. A much cleaner approach is 
to fully exploit the symmetries of the problem, which is 
what we shall do. 

We explained in appendix [X] that generalized spin 
wavefunctions are classified by irreducible representa- 
tions of the Lie algebra of SU(n) given by forming tensor 
products of N fundamental multiplets. Specifically, we 
can decompose the vector space of generalized spin wave- 
functions into sub-spaces labeled by Young tableaux of 
shape A and in such a decomposition the multiplet la- 
beled by A occurs / A times. We can thus write a basis of 



generalized spin wavefunctions X gsw (TV, A, A z , r) labeled 
in the same way as states within SU(n) multiplets, where 
TV is the number of particles, A is the symmetry type, 
and where r runs from 1 to f x and distinguishes between 
SU(n) multiplets occurring in the tensor product with the 
same A. We shall use the convention that wavefunctions 
with different r values are mutually orthogonal. Recall 
from section IIII CI (and see below in section IB ip that 
A z describes how many particles in the generalized spin 
wavefunction are of each type a, f3, 7, . . .. Note that it 
may or may not be possible to construct a wavefunction 
in representation A given a particular TV and A z as we 
will describe further in section [B~T1 below. 

Given the above decomposition of the generalized spin 
wavefunction into representations, we can analogously 
decompose the combined wavefunction ip based on the 
same irreducible representation of SU (n) . This follows 
because the spatial function is independent of the the 
component values. Thus, without loss of generality, we 
may divide the space of possible TV-particle wavefunc- 
tions ip into subspaces indexed by A and A z . In or- 
der to decompose our wavefunctions into a basis of ir- 
reducible representations whilst remaining of the general 
form given in Eq. IB2l wc propose the following alternative 
general form for a basis for ip, which loses no generality: 



f 



CrA $ S patialY gsw (TV, A, A Z ,T*) , (B3) 



where c r are arbitrary coefficients. For bosons wc simply 
use a symmetrization operator S instead of an antisym- 
metrization operator A. 

In this expression we are free to choose the most gen- 
eral possible $ S patiai- However, here not all different 
choices of ^spatial will correspond to different wavefunc- 
tions ip. In the section IB II we will first construct the 
generalized spin wavefunctions X gsvl (TV, A, A z , r); then in 
section IB 21 we will (use group theory to) determine how 
to construct a complete (but not over-complete) basis for 
the spatial wavefunctions $ S patiai- 



1. Introduction to Generalized Spin Wavefunctions 

The generalized spin wavefunctions X gsw (TV, A, A z ,r) 
describing a system with n components correspond to 
particular irreducible representations of SU(n): those 
given by forming tensor products of N fundamental mul- 
tiplets. Such representations are in correspondence with 
the standard tableaux of shape A and hence with a subset 
of the symmetric group. We must specify the tableaux 
shape A, where A is a partition of the integer N into 
at most n parts, the standard tableaux index r, and the 
component composition of the function, which is denoted 
by X z . In order to enforce that X gsw (TV, A, A z , r) are mu- 
tually orthogonal in A, A z and r, we shall construct our 
SU(n) representations based on the orthogonal represen- 
tation of the symmetric group. 
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The generalized spin wavefunctions can be written as 
linear combinations of primitive generalized spin wave- 
functions: these are defined by specifying a list of the 
the values that each component can take; a, /?, 7 and so 
on. For example = |a7/3/3a) is a primitive generalized 
spin wavefunction. One can order these wavefunctions 
according to a lexicographic ordering scheme, so that the 
first primitive generalized spin wavefunction would be 

tfi (TV, X z ) = \aaa...j3/3...CO ■ 

The index X z represents the list {TV aj TV^, TV^} of TV 
which corresponds to the number of times each type 
of component value occurs in the primitive generalized 
spin wavefunction. In general we have wavefunctions $j , 
where j runs from 1 to the number of primitive gener- 
alized spin wavefunctions d(X z ) of all possible orderings 
of the components, and indexes the position of $j within 
an ordered list of those wavefunctions. 

The explicit forms of the generalized spin wavefunc- 
tions can be deduced as follows: first we write down the 
orthogonal representation matrices U(iTi) x correspond- 
ing to the subset of standard tableaux of shape A; from 
this orthogonal representation we then construct a corre- 
sponding group algebra using the matric basis, Eq. IA1I 
As the elements of the group algebra are expressed as 
a linear combination of permutation operators, one can 
apply the matric unit operator e A s to any primitive gen- 
eralized spin wavefunction to give either zero, or a lin- 
ear combination of permuted primitive generalized spin 
wavefunctions. A key result of this construction is that 
there is an irreducible representation of the symmetric 
group associated with each A and that representations 
associated with different A are inequivalent^. Thus the 
generalized spin wavefunctions constructed in this way 
are in one-to-one correspondence with the irreducible 
representations of SU(n) and with the SU(n) multiplcts 
occurring in the tensor product of TV fundamental multi- 
plcts. In other words, the operator e A s acts as a projector 
on the primitive generalized spin wavefunctions to give a 
particular representation of SU(n). 

A further key result is that the wavefunctions 
e^ s -dj(N, X z ) are either zero or, independently of the 
choice of r and s, they belong to the same symmetry 
type A. Due to the orthogonality of the matric unit oper- 
ators, Eq. IA21 two such wavefunctions e^ s dj(N, X z ) and 
e^, s t9j(TV, X z ) are orthogonal if r^r'. The operators e A s 
act as shift operators in the sense that 

e A s (e A s ^(TV,A,)) = e A s ^-(TV, A z ), 

and so operating on the first projected function e A s #j 
yields a set of orthogonal functions belonging to the same 
symmetry type , and further, by acting with all such 
possible shift operators on all possible -dj, we obtain a 
complete basis of states having the same A and X z . 

For the SU(2) case it is possible to prove further that 
S 2 e^Jj(N, S z ) = S{S+l)e^j(N, S z ) with A the repre- 
sentation corresponding to spin S (the proof is given in 



Ref. In other words, for the two-component case, the 
generalized spin wavefunctions are precisely spin eigen- 
f unctions. 

Thus we may decompose the complete space of gener- 
alized spin wavefunctions X gsw in terms of the multiplet 
label A and r, and the multiplet state label X z . Most 
generally, we write 

X gsw (TV, A, A„r) = £ £ b 3j e A s ^(TV, A z ), (B4) 

8=1 j=l 

where 6 SJ are arbitrary coefficients. Since these wave- 
functions are constructed in correspondence with the 
multiplets of SU(n), the number of linearly independent 
values of b s j is precisely the multiplicity, M(X, X z ), of the 
SU(n) multiplet state (this can be explicitly checked). It 
can be shown that these wavefunctions arc mutually or- 
thogonal in A, X z and r due to the orthogonalit y re lations 
satisfied by the matric unit operators (see Ref. |64| for the 
proof in the SU(2) case). 

The mapping to representations of the symmetric 
group manifests itself in the permutation symmetry of 
the wavefunction, for example TV gsw (TV, A, X z , 1), inde- 
pendent of A z , must be symmetric in labels 1 to TVi, in 
labels TVi + I to TV2 and so on where TV^ are the integers 
describing the partition A. We shall describe the explicit 
forms of the generalized spin wavefunctions in more de- 
tail in appendix [Cj 

2. Combination of Generalized Spin Wavefunctions 
and and Spatial Wavefunctions 

In this section of the appendix we shall now describe 
the procedure for sewing together the spatial and general- 
ized spin wavefunctions to give a combined wavefunction 
with the correct symmetry property. We shall then focus 
on the explicit form of the spatial part of the wavefunc- 
tion. Wc aim to determine the dimension of the vector 
space of spatial wavefunctions. 

So far we have shown that the most general expression 
for the basis of combined generalized spin wavefunctions 
is given in Eq. IB3I It is insightful to split up the an- 
tisymmctrization operator into parts acting separately 
on either the spatial or spin functions — a permutation 
TTi is equivalent to permutations irf acting on the spa- 
tial wavefunction and 7rf acting on the generalized spin 
wavefunction. Thus Eq. IB3I can be written in the equiv- 
alent form: 

^(TV,A,A Z ) 

f X N\ 

= Yj <vE s S n ( 7r 7r f *s P atiai7rf A" gsw (TV, A, X z ,r) . 

r=l i=l 

Now, the most general expression for the generalized spin 
wavefunctions is given in Eq. IB4I An important result, 
which follows from the definition of the matric unit op- 
erators, is that the action of an arbitrary permutation 
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operator irf on a generalized spin wavefunction is: 



Trf X gsw (AT, X,\ z ,r)=J2 ^g- (N, A, X z ,s) U(tt^ 



Using this result, we can write Eq. IB3I without the anti- 
symmetrization operator in the form: 

ip (N, A, X z ) = c r X ESW (N, A, X z , s) $ sr , 



where we have introduced the spatial wavefunction $ sr , 
which is defined as: 

= ( Jn ) Z ^(^Ji-sgn fa)*? $ spatial . 

Finally, using the definition of the contragradient matric 
unit operators, we can write the spatial wavefunction in 
the form 



m 



1/2 



spatial ■ 



In the boson case, where we require the wavefunction to 
be fully symmetric (thus removing the sgn (jTi) factor), it 
clearly follows that A is replaced by A. Note that $ spatial 
is still given by Eq. IB11 

Now, in order to simplify our presentation, we shall de- 
fine the rth part of spatial wavefunction to be the spatial 
wavefunction associated with the rth multi-component 
function, in other words: 



f 

$ (N, A, r) = ^ c s e^$ spatial , 

s=l 



(B5) 



or with A replaced with A for the boson case. These 
spatial wavefunctions are mutually orthogonal in A and 
in r due to the othogonality relations satisfied by the 
matric unit operators. Note that the explicit form of 
$ (N, A, r) will be: 

JV orbitals 

$(N,X,r) = J2 (»*))■ (B6) 

»ll...tjv = l 

where now the coefficients a' depend on the combinations 
of the a coefficients (from Eq. IBI[) and the c coefficients 
(introduced in Eq. IB3[) and in particular they must con- 
tain linear dependences. 

Thus, the most general form of the basis wavefunctions 

is 

/ A 

$ (N, A, A 2 ) = J2 X s™ (N, A, A 2 , r ) $ (N, A, r ) , (B7) 

r=l 

Since both X gsw (N, A, X z , r) and ip(N 7 X,r) are mutu- 
ally orthogonal in A, A 2 and in r, the wavefunction 



ip (N, A, X z ) is itself othogonal to wavefunctions with dif- 
ferent A and X z values. We now propose to simplify this 
result using the following theorem. 

Theorem 1: The wavefunction ip(N,X,X z ) 
in Eq. \B7\ can be written in the form: 

A&(N,X,l)X gsw (N,X,X z ,l), for fermions, or 
<S$ (N, A, 1) X gsw (N, A, A 2 , 1) for bosons. 

Proof: Let us demonstrate the proof in the fcrmion case 
only, since the boson case follows along very similar lines. 
The proof comes from the following two results: firstly 
from Eq. IB41 using the orthogonality condition for the 
matric unit operators, we have: 

*g SW (N, A, A 2 , r ) = £ b Sij e^j (N, X z ) , 

s=l,j=l 

= eri E " Ki ZisMN, A 2 ) - e^X ssw (N, A, A 2 , 1). 

Thus all coefficients in X gsw (N, A, A 2 , r) arc linearly de- 
pendent on the set of coefficients in X gsw (N 7 A, A 2 , 1). 
Similarly, the second result is that 

^ ' C s e rs $ S p a tial = C r i ^ C s e ls $ S p a ti a l, 
s=l s=l 

or alternatively 

$(N,X,r) = e$ 1 $(N,X,l), 

and so the coefficients in 4> (N, A, r) are linearly depen- 
dent upon the coefficients in $ (AT, A, I). 

Now substituting these two results into Eq. IB7I we 
have: 

/ A 

V (N, A, A 2 ) = J2 (N, A, I) e^X gsw (N, A, A 2 , 1). 

r=l 

Then using the definitions of the matric unit operators, 
this simplifies to 

Nl 

1& (N, A, A 2 ) cx J2 sgn (tt 4 )^ ($ (N, A, I) X gsw (N, A, A 2 , 1)) 



r=l 

and since the U(tti) x matrices are orthogonal matrices, 
this further simplifies to 

^(N, A, A 2 ) cx i$ (N, A, I) X gsw (N, A, A 2 , 1), (B8) 

and thus we have proved theorem 1. For the boson case, 
we replace the contragradient representations A by A, and 
consequently we will arrive at a symmctrizer S instead 
of the antisymmetrizer in the final expression. 

Corollary: We have expressed all of the linearly inde- 
pendent terms in the basis wavefunctions ip(N, A, A 2 ) in 
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terms of the product $ (N, A, 1) X gsw (N, A, X z , 1) — the 
entire basis wavcfunction can then be recovered from this 
single term by simply antisymmetrizing or symmetrizing 
in all the particle co-ordinates, and further, this anti- 
symmctrization procedure does not reduce the number of 
linearly independent terms. In other words, the dimen- 
sion of the space of multi-component basis wavefunctions 
ip (N, A, X z ) is given by the dimension of the space of gen- 
eralized spin wavefunctions X gsvl (N, A, A z , 1) multiplied 
by the dimension of the space of spatial wavefunctions 
$ (N, A, 1). Our task now is to determine the dimensions 
of these vector spaces. We shall discuss the generalized 
spin wavefunctions X gsw (iV, A, X z , 1) for various cases in 
appendix [C] and we shall discuss the spatial wavefunc- 
tions $ (N, A, 1) in appendix [DJ 



Appendix C: Vector Space of Generalized Spin 
Wavefunctions 

In this appendix we shall explain how to derive the 
space of generalized spin wavefunctions X gsw (N, A, X z ,r) 
in a variety of cases. First we shall explain the simplifica- 
tions that arise in the description of the two-component 
case (spin-1/2 eigenfunctions) ; we shall then explain the 
most general n-component case; in the final two sections 
we shall consider two special cases — the description of 

3- component systems with spin-1 and the description of 

4- component systems with spin-1/2 combined with a two- 
valley or bi-layer degree of freedom. In each case we shall 
describe how the procedure described in appendix [B] is 
applied to generate the vector space of generalized spin 
wavefunctions. 



1. Two-Component Systems 

An important observation is that any two-component 
multi-particle system can be mapped onto a system with 
a spin-1/2 degree of freedom regardless of whether it con- 
tains fermionic or bosonic particles. The procedure of de- 
composing into SU(2) representations is then identical to 
the problem of the construction of spin eigenfunctions — 
a well studied problem in theoretical chemistry^. Note 
that in general, a wavcfunction need only be an eigen- 
function of spin if the system is invariant under spin ro- 
tations. Nonetheless, since a decomposition in terms of 
spin eigenfunctions is akin to a decomposition in terms 
of SU(2), the spin eigenfunctions will always provide a 
complete basis in which to describe any two-component 
system. 

Spin eigenfunctions for a multi-particle system are de- 
fined to be eigenfunctions of the total spin operator S 2 
and total z-component of spin operator S z . Such func- 
tions can be expressed as linear combinations of primi- 
tive spin wavefunctions : to repeat the definition given in 
section fill Al a primitive spin wavefunction $j of a many 
particle system is an eigenfunction of the S z operator of 



every particle in the system. (They are effectively two- 
component generalized spin wavefunctions where a = ||) 
and P = ||). ) 

Spin eigenfunctions X(N,S,S z ,r) are constructed 
from linear combinations of these primitive spin wave- 
functions. For a system of N particles with spin 
eigenvalue S, there is a one-to-one correspondence be- 
tween the spin eigenfunctions and the subset of or- 
thogonal irreducible representations of the symmet- 
ric group denoted by Young tableaux of the shape 
A = [i N + S, \N - S] , see Ref. M 

The following results apply specifically to the two- 
component case. For any primitive spin wavefunction 
dj it can be shown that 

e^ r i9j oc e*i#i, or 0. 

A corollary is that one can generate the entire space 
of spin eigenfunctions from only the first primitive spin 
wavcfunction d\. The key expression for the construction 
of spin eigenfunctions is then 

X{N,S,S z ,r)o^e\^. (CI) 

Let us demonstrate these ideas with a simple example: 
the 3 particle state with eigenvalues 5=1/2 and S z = 
1/2. The associated Young tableaux are of the shape 
A = [2, 1] : 



There are two possible arrangements of standard tableau: 



1 


3] 


1 


2| 






3^ 





and there are correspondingly two eigenfunctions with 
the same spin eigenvalue, each of which is generated by 
one of the group algebra basis operators 

A(3,i,i,l) =4 2 1 ' 1 V 1 (3,1/2) 

« 7H (2 im> - mt) - im» 
x (3,i, 1,2) = gp' 1 ]^ (3, i/2) ex ^ (int> - \m) ■ 

Upon inspection we see that X (3, |, \ , l) is symmetric 
under exchange of indices 1 and 2 with no particular 
symmetry conditions for index 3, also X [3, 5,5,2) is 
antisymmetric under exchange of the same indices. 

Recall the most general expression for the wavefunc- 
tion given in Eq. IB8I In the two-component case, since 
X(N,X,X Z ,1) = X(N, S,S„1) = e^MN,S z ), we need 
only specify the first primitive spin wavcfunction in or- 
der to describe the full wavefunction. Thus in the two- 
component case the wavcfunction (for fermions) can be 
succinctly expressed as: 

i/>(N, S, S z ) oc M {N, S, 1) i?i(JV, S z ). (C2) 
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2. Non Spin-Rotationally Invariant n-Component 
Systems 

In the most general n-component case, the system may 
not be invariant under spin rotations, or, if there is no 
spin degree of freedom present, rotations within the gen- 
eralized spin space. The most general set of generalized 
spin wavefunctions in this case are simply those specified 
by the decomposition into tensor products of N funda- 
mental SU(n) multiplets. The irreducible representations 
occurring in the decomposition are labeled by A and r, 
with no additional constraints — these are the functions 
described by Eq. IB4I We note that n-component sys- 
tems of this type have previously been conceived in Rcf. 
l67l . however the explicit wavefunctions have not been cal- 
culated. 

Since primitive generalized spin wavefunctions at dif- 
ferent positions in the ordered set are related by 
permutations, it is possible to write Eq. IB4I (i.e., 
-^gsw (N, A, A 2 , r)) in the form: 

/\d(A„) AI(\,\ Z ) 

S—1J—1 p—1 

(C3) 

where p runs from 1 up to the multiplicity, M(A, A 2 ), 
of the SU(n) multiplet state; recall that this is given by 
the number of semi-standard tableaux of shape A and 
component content A 2 . Note that since p can be greater 
than 1 but is always less than or equal to d(X z ) (that 
is, the total number of primitive generalized spin wave- 
functions including all possible ordcrings of components), 
we include only a certain subset of the lexicographically 
ordered primitive generalized spin wavefunctions in the 
basis — these correspond precisely to the semi-standard 
tableaux shapes. For example, in the A = [3, 1] state of 
SU(3) the semi-standard tableaux are 
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The corresponding primitive generalized spin wavefunc- 
tions are given by interpreting the placement of the num- 
bers 1 to n in the semi-standard tableaux as the place- 
ment of the components a, /?, 7, . . . in the primitive gen- 
eralized spin wavefunctions. The corresponding set of 
primitive generalized spin wavefunctions that can occur 
in this state are then: 

\aaaj3) , \aaaj) , \aa/3j3) , , |aa 7l 9) , 

|aa 77 ) , \a000) , |a/3/3 7 ) , |a/3 7 /3> , |a/3 77 ) , 
|a 77 /3) , |a 777 ) , \P00y) , |/3/3 77 ) , |/? 777 ) • 

The set of & P {N,X Z ) for TV up to 5 and for SU(3) and 
SU(4) arc listed in tabic HVl The results presented in ta- 



bic IIVI have been checked by explicitly applying the ma- 
tric unit operators to all possible primitive generalized 
spin wavefunctions and then determining the linearly in- 
dependent combinations algebraically. The results of this 
brute-force approach for the numbers of linearly inde- 
pendent state precisely match with the multiplicities and 
multiplet dimensions deduced by the decomposition into 
irreducible representations of the Lie algebra of SU(n) 
according to appendix [X] 

Recall that in the tensor product of N fundamental 
SU(n) multiplets there are / A irreducible representations 
associated with each symmetry type A, for example with 
A = [2, 1] we have f x = 2, and thus there are / A gener- 
alized spin wavefunctions of the form X gsw (N, A, X z ,r). 
When we construct the fully symmetric/antisymmetric 
basis wavefunction by sewing together spin and spatial 
parts (Eq. IB3[) . any orthogonal generalized spin wave- 
functions of the same symmetry type and generated by 
the same root primitive generalized spin wavefunction 
will appear in the same basis wavefunction. Thus, when 
counting the number of generalized spin wavefunctions 
we include this degeneracy (see table IIV[) . but when 
counting the size of the vector space of generalized spin 
wavefunctions wc do not . 

Recalling Eq. IB81 we aim to determine the most gen- 
eral form for X gsvl (N, A, A 2 , 1). We can use Eq. IC3I to 
write the basis of X gsw (N, A, A 2 , 1) in terms of ap- 
plied to the set of primitive generalized spin wavefunc- 
tions given in tabic IIVI Thus the minimal information to 
reconstruct the wavefunction is contained within that set 
of primitive generalized spin wavefunctions. In the multi- 
component case the wavefunction (for fermions) can then 
be succinctly expressed as: 

iftN, A, X z ,p) oc M (N, A, 1) P (N, A 2 ), (C4) 

where the set of functions $j,(iV, A 2 ) are taken from table 

El 

3. Spin-Rotationally Invariant n-Component 
Systems 

If it transpires that the system which we are describing 
is spin rotationally invariant, then it is convenient to con- 
struct a basis of spin cigenfunctions; the Hamiltonian will 
be diagonal in such a basis. We propose to construct ap- 
propriate spin eigenfunctions by forming particular linear 
combinations of the generalized spin wavefunctions given 
in Eq. IB4I The resulting wavefunctions will be labeled 
by spin eigenvalues S and S z . Let us begin by examining 
the origin of the basis in more detail. 

A single particle which is symmetric under three- 
dimensional spin rotation is invariant under the three- 
dimensional rotation group 0(3) A These irreducible 
representations of the rotation group are characterized 
by the standard set of spin angular momentum quan- 
tum numbers, j and j z . There exists a general map- 
ping between such irreducible representations of the ro- 
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tation group and irreducible representations of the group 
SU(2j + an d in particular there is a one-to-one 

mapping between a system of particles which individ- 
ually have spin angular momentum j = 1/2 and irre- 
ducible representations of SU(2) (since SU(2) is isomor- 
phic to 0(3) ). In addition there exists a further one- 
to-one mapping between irreducible representations of 
SU(2j + 1) and irreducible representations of the sym- 
metric group, specifically the irreducible representations 
of SU(2j+l) visualized by Young tableaux containing N 
boxes in no more than 2j+l rows: A = [N±, N2, A^'+i] 
with iVi > N 2 > ... > N 2] +i > 0. 

Our general approach has been to decompose general- 
ized spin wavefunctions into a basis built from irreducible 
representations of the group SU(n). Due to the isomor- 
phism with the rotation group for the case of n = 2, we 
find that our general approach in this case corresponds 
precisely to the construction of spin cigenfunctions. For 
n > 2 however the mapping becomes a little more com- 
plicated. In order to write our generalized spin wave- 
functions in a basis of spin cigenfunctions what we are 
really doing is decomposing representations of SU(n) into 
representations of 0(3). The procedure for this decom- 
position is explained in Ref. |4g, chapter 11: the principle 
feature of this result is that this decomposition is not a 
one-to-one mapping for a general SU(n,) decomposition, 
in other words, due to this mapping, it follows that there 
are multiple spin eigenvalues corresponding to each A. 
We must therefore label the basis by both A and by the 
spin eigenvalues S and S z , that is, our set of spin eigen- 
functions are X (N, A, S, S z ,r) i.e., we have now replaced 
X z with S and S z . These functions remain orthogonal in 
A and in r. 

We shall now explain our procedure for how to con- 
struct spin eigenfunctions that describe spin-1 particles 
(the general case can be deduced along similar lines, and 
checked using the information contained in Ref. l48h . The 
explicit forms of the spin eigenfunctions can be deter- 
mined as follows: first associate a spin value with our 
notation such that, for each particle, a represents the 
s, = 1 state, /3 represents the s z = and 7 represents 
the s z = — 1 state; we note that there could be several A 2 
values that could correspond to the same S z value, and so 
we form a linear combination of these states, for example 
A 2 = {1, 1, 1} and A 2 = {0, 3, 0} have the same S z value 
0, and so we must work with the linear combination 

X (N, A, S, S z = 0, r) 
= AX gsw (N, A, {1, 1, 1} , r) + BX gsw (N, A, {0, 3, 0}],r) , 

where A and B are arbitrary coefficients. In general we 
would form a linear combination of all possible A z for a 
given S z . Then we apply the S 2 operator and solve the 
self-consistency condition: 

S 2 X (N, A, S, S z , r) = X (N, A, S, S z , r) , (C5) 

for the coefficients A, B, 



We shall demonstrate our technique with some simple 
examples. Let us first consider N — 3, A = [3], S z = 1 
and r = 1. In this case the valid primitive generalized 
spin wavefunctions are: 

\aaj) , \a/3/3) . 

Applying S 2 to this vector space gives: 

a 2(3 [3] ( \aary) \ _ ( 4 4 \ 43] / \aa-y) \ 

Using this information we now apply Eq. IC5I and we 
arrive at 

S 2 e [ ?l(A\aaj) +B\af3f3)) 

= efl((4A + AB) \aa<y) + (44 + 10 B) 

= S{S + l)ef }(A \aaj) + B | a/3/3)). 

Hence we must solve the eigenvalue equation 

(Notice that the matrix appearing in the eigenvalue equa- 
tion is the transpose of the matrix written in the first 
step.) The solutions for X (N, A, S, S z , 1) in this case 
correspond to S = 1,3, and the associated eigenvectors 
are 

eg (|aa7> + 2 \affi) , 
(|a/3/3) - 2 |aa 7 » , 

which are now precisely spin eigenfunctions. Given these 
results once can obtain the set of orthogonal spin eigen- 
functions X (AT, A, S, S z ,r) by applying the ef\ operator. 

The solutions obtained here (S = 1,3) are precisely 
the spin eigenvalues corresponding to the decomposition 
of SU(3) into 0(3) as described in Ref. 0. We have de- 
liberately constructed the eigenfunctions in such a way 
that they can be specified by a linear combination of 
primitive generalized spin wavefunctions X\ c (for exam- 
ple X\ c = |a/3/3) — 2 I aery)) with a projection operator 
applied to them: 

X(N,X,S,S z ,l) = e^X lc . 

The projection operator does not depend on the set of 
spin eigenvalues, only on the symmetry type. Appropri- 
ate forms for X\ c are given in table [V] Note that table 
fVl only lists cases where S — S z : spin eigenfunctions for 
S z < S can be deduced by explicitly applying a spin 
lowering operator S~ . We have derived the results pre- 
sented in table [V] by precisely the procedure outlined in 
this appendix. 

Recall Eq. IB8I In the case described in the section, the 
minimal information required to construct the full wave- 
function is just X (N, A, S, S z , 1). Since these functions 
can be expressed in the form e^ATic, we have deduced 
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that minimal information to reconstruct the wavcfunc- 
tion is in fact contained within these linear combinations 
X\ c . Thus the full wavefunction (for fermions) can be 
expressed as: 

1>(N,\,S,S,)<xM (N,X,l)X lc) 

Let us give one further example briefly. Consider the 
N = 4, A = [3, 1], S z = 1 state. In this case we must in- 
clude the following set of primitive generalized spin wave- 
functions: 

\aah) , \aar/P) , \a0ffl 
Performing the same analysis we find: 

/ \aaPi) \ I 10 2 \ / \aafii) \ 

S 2 e [ ^ ] \acrff}) = 2 4 2 ef^ \aorfP) 

\ \am) ) V 4 4 6 J \ \ a ppp) j 

In this case we find that the possible S values are 3,2 
and 1, precisely in correspondence with the underlying 
decomposition of SU(3) representations into 0(3) repre- 
sentations. The resulting spin eigenfunctions obtained 
for this case are given in table [Vj 

More generally for a system of spin j particles (an ex- 
ample of an n = 2j + 1 component system) we can con- 
struct a basis of spin eigenfunctions by decomposing rep- 
resentations of SU(2j+l) into representations of 0(3), us- 
ing a similar enumeration procedure to the 3-component 
case discussed above. Tables of such decompositions are 
given in Ref. for states up to j = 7/2. 



two valleys labeled by 1 and 2 or A and B. The same 
situation could arise in a bilayer system, where there 
could be two layers (1 and 2) of spin- 1/2 particles. Al- 
though we are dealing with four-component wavefunc- 
tions here, the spin component itself only corresponds to 
spin-1/2 (as opposed to the four-component spin-3/2 sys- 
tem say). If the system is spin-rotationally invariant, but 
not invariant under valley symmetry, then it is fruitful to 
construct a basis in terms of spin eigenfunctions for the 
two types of valley electrons. As before, we propose to 
construct appropriate spin eigenfunctions by writing lin- 
ear combinations of the generalized spin wavefunctions 
which have the same S z eigenvalue and then solving a 
self-consistency equation as in Eq. IC5I 



A four-component system has components labeled by 
a, p, 7 and S. In order to write down spin eigenfunctions 
we specify the mapping: in valley 1, a corresponds to the 
s z = 1/2 state and j3 corresponds to the s z = —1/2 state, 
and in valley 2, 7 corresponds to the s z = +1/2 state and 
8 corresponds to the s z = —1/2 state. We can then pro- 
ceed to construct spin eigenfunctions along the same lines 
as in appendix IC 31 first we enumerate all possible lin- 
early independent four-component generalized spin wave- 
functions; then we apply the S 2 operator to each; finally 
we determine linear combinations of the four-component 
generalized spin wavefunctions which result in spin eigen- 
functions. A more rigorous mathematical underpinning 
of our method in terms of representation theory is pre- 
sented in appendix IE 31 



4. Partially Spin-Rotationally Invariant 
n-Component Systems (Spin+ Valley or Spin+Layer) 

A system such as graphene is a four-component sys- 
tem: there are two types of electrons corresponding to 



As an example, consider the 3-particle state with sym- 
metry type [2,1] and S z = 1/2. Applying the S 2 operator 
to the vector space of generalized spin wavefunctions we 
find: 
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Using this information we implement Eq. IC5I the solu- 
tions for the spin eigenfunctions are given by determining 
the eigenvalues and eigenvectors of the matrix appearing 
here. We find in this example that 5 = 3/2 occurs twice, 
and 5=1/2 occurs 6 times. Due to this multiplicity we 



are forced to introduce an extra index, k, to distinguish 
between these linearly independent and orthogonal states 
with the same S eigenvalue. Our spin eigenfunctions are 
now written as X (TV, A, 5, k, S z ,r). As in appendix IC 31 
the solutions for the spin eigenfunctions (given by the 
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eigenvectors of the above matrix) can be written as a 
projection operator applied to a given linear combina- 
tion of primitive generalized spin wavef unctions: 

X(N,\,S,k,S z ,l) = e$ 1 X lc . 

The appropriate linear combinations X\ c of primitive gen- 
eralized spin wavefunctions arc given in table IVII Note 
that table IVII only lists cases where S = S z : spin eigen- 
functions for S z < S can be deduced by explicitly ap- 
plying a spin lowering operator S~ . We have derived 
the results presented in tabic PVTl bv precisely the proce- 
dure outlined in this appendix. Also, following the same 
reasoning given in appendix IC 3[ the minimal informa- 
tion to construct the wavefunction is contained within 
the expressions X\ c , and thus the full wavefunction can 
be re-constructued from 

V>(TV, A, S, k, S z ) oc i$ (TV, A, 1) X lc . 



Appendix D: Spatial Wavefunctions for Quantum 
Hall Systems 

In this appendix we will apply the general techniques 
for the construction of spatial wavefunctions derived in 
appendix [B] to the specific case of translationally invari- 
ant LLL wavefunctions of fixed degree L. Specifically 
out task is to determine the functions $ (TV, A, 1) which 
will provide the minimal information required to recon- 
struct the complete wavefunction. The enumeration of 
these functions will lead us to a derivation of the dimen- 
sion of the space of spatial wavefunctions. In appendix 
ID II we shall describe how to construct appropriate LLL 
wavefunctions in general, and we shall explain how the 
vector space of LLL spatial wavefunctions can be enu- 
merated algebraically using a set of projection operators. 
In appendix ID 21 we shall explain a more general proce- 
dure to calculate the dimensions of the vector space of 
LLL spatial wavefunctions. 



1. Procedure for Construction of Spatial 
Wavefunctions in the LLL 

Spatial wavefunctions describing multi-particle LLL 
states are built out of sums of products of single par- 
ticle LLL wavefunctions and as we argued previously, at 
the beginning of section [Til it follows that the spatial 
wavefunctions must be (translationally invariant) ana- 
lytic polynomials of a fixed degree in the complex rela- 
tive co-ordinates Zi. For a state with fixed relative an- 
gular momentum eigenvalue L the most general form of 
spatial wavefunction is given by: 

JV 

^spatial = ^ a il ... iN zl 1 ... z l £ with y^ i m = L. 

ii,...ijv m— 1 



with arbitrary coefficients ai 1 ...i N (c.f. Eq. IB1[) . For 
example, for TV = 3 and L = 2 the most general spatial 
wavefunction is: 

^spatial = 0-200^1 + a 020^2 

+aoo2^| + ai W zih + a nz 2 z 3 + a wl ziz 3 - 

Using this definition wc could in principle now use Eq. 
IB5I to produce an expression for <3> (TV, A, 1) in the form 
given in Eq. IB6I In practice, we can employ a further 
short cut to this procedure that reproduces exactly the 
same result. Recalling the proof of theorem 1, we showed 
that: 

$(TV,A,r) = e^$(TV,A,l). 
This result implies the following idempotency relation: 

$(TV,A ! l)=ei 1 $(TV ! A,l). 

Thus we see that if we start by writing $ (TV, S, 1) as the 
most general form of spatial wavefunction ^spatial given 
above, then the action of applying e\ x will automatically 
project out the most general form of the first spatial 
wavefunction, given in Eq. IB6I We can interpret this 
process of applying the symmetric group algebra basis 
operator as a process of introducing many linear depen- 
dences between the original set of linearly independent 
coefficients a,i^...i N - The result will be a polynomial of 
a particular symmetry imposed by the representation of 
the symmetric group A, which may still contain some lin- 
early independent coefficients a\ ijv — the number of 
which is the dimension of the space of spatial wavefunc- 
tions. 

Wc have now deduced that the vector space of LLL spa- 
tial wavefunctions can be projected out from the space 
of most general LLL wavefunctions. In principle this 
method can be employed to determine the dimensions of 
the basis (listed in tabic UJ and the explicit forms of the 
basis functions. However, due to the algebraic complex- 
ity of constructing the most general form of the spatial 
wavefunctions, and the fact that the number of terms 
in e^i grows as TV! with the number of particles TV, this 
brute force approach is not viable in general. In the next 
section we shall describe a method to calculate the di- 
mensions of the space of LLL spatial wavefunctions in 
general. 

2. Calculating the Dimension of the Vector Space 
of Spatial Wavefunctions 

In this section we shall explain the method used to cal- 
culate the values listed in table UJ Using these results we 
shall then re-interpret the form of the spatial wavefunc- 
tions $ (TV, A, 1) and in doing so we shall complete our 
explanation of the results presented in table [UJ 

The spatial wavefunciton $ (TV, A, 1) takes the form of a 
translationally invariant, analytic polynomial in the co- 
ordinates zi, satisfying a particular set of permutation 
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symmetries imposed by the underlying symmetric group 
representation. In general it contains many linearly in- 
dependent coefficients, the number of which depends on 
the order of the polynomial. This is called the dimen- 
sion of the polynomial space. In order to determine how 
many independent parameters occur at a given order we 
construct generating functions for the polynomial dimen- 
sion, and to do this we first need to construct a linearly 
independent basis in which to describe the polynomials. 

Recall from section [TT] that any antisymmetric polyno- 
mial in N variables can always be written as a symmet- 
ric polynomial by factoring out a Jastrow factor, equation 
([7]) . There are many ways to construct a basis of symmet- 
ric polynomials; the one we utilize is the basis of elemen- 
tary symmetric polynomials, which we defined in equa- 
tion l|8|). If we impose the condition that the polynomi- 
als must be translationally invariant then we must write 
elementary symmetric polynomials in terms of relative 
co-ordinates i,-, which we defined in equation (jH]). Recall 
also that translational invariance results in equation (TT0|) . 
These modified elementary symmetric polynomials form 
a basis of translationally invariant symmetric polynomi- 
als, in other words we can form any translationally in- 
variant symmetric polynomial from a sum of products of 
the modified elementary symmetric polynomials. (More 
precisely, these polynomials form what is called a ring of 
translationally invariant symmetric polynomials^.) 

The polynomials such as $ (N, A, 1) are not fully sym- 
metric, but instead they have symmetries in subsets of 
particle indices as dictated by the shape A of their asso- 
ciated Young tableaux. As a starting point we shall try 
to use the basis of elementary symmetric polynomials to 
span the space of polynomials with subsets of symme- 
tries. We require in addition that the polynomials are 
translationally invariant and this places a restriction on 
the basis, 



ei,i...Ar 



0. 



(Dl) 



The dimension c? S ym -^0 of the space of translation- 
ally invariant symmetric polynomials in TV variables and 
of degree L, is defined in terms of the generating function 
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9=0 



We can write down a generating function for the space 
of polynomials with subsets of symmetries by multiply- 
ing the generating functions of the subsets (for proof see 
appendix IE 4|) . 

Our task now is to construct a basis of partially sym- 
metric polynomials which reflects the symmetric group 
representations. It happens that the two component 



case is the simplest to explain, and so we shall discuss 
that case first before generalizing our result to the multi- 
component case. 



a. Two-Component Case 

The basis of two-component wavefunctions is identical 
to the basis of spin eigenfunctions. We shall thus consider 
the spatial funcitons associated with an N particle state 
of spin S. 

For bosons the particular case of interest is a transla- 
tionally invariant polynomial which is symmetric in ar- 
guments 1 to + S and in the arguments + S + 1 
to N. Due to equation (|Dip we must eliminate elemen- 
tary symmetric polynomials of degree 1 from the gen- 
erating function of one subset. The generating function 
for translationally invariant, symmetric polynomials with 
two symmetry subsets is then 
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For example, we might consider polynomials of order 
L = 3 in N = 5 indices which are symmetric in indices 
1,2 and 3, and separately symmetric in indices 4 and 5. 
In terms of elementary symmetric polynomials we can 
construct: 



61,12362,123 

e3,123 



ei, 12362. 45 
p 3 

e l,123 



This number, 4, is precisely given by appropriately dif- 
ferentiating the generating function given in Eq. ID31 for 
the case of S = 1/2. 

For fermions we require that the polynomial is sep- 
arately antisymmetric in the two symmetry subsets in- 
stead of separately symmetric in the two subsets, and so 
the space of symmetric polynomials only comes into play 
at polynomial orders higher than the order of the Jastrow 
factor. To start the generating function at a particular 
order, J, we must multiply by q J (for a proof see ap- 
pendix IE 4p . In order to take into account the Jastrow 
factors associated with each of the symmetry subsets we 
must have 
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The generating function for translationally invariant, 
antisymmetric polynomials with two symmetry subsets 
is then 



N 



— S 



N 



-S 



^fcrmion 
J N,S 



q J - J - 1 — q" 

m=2 H n=l H 
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More precisely these generating functions actually ac- 
count for polynomials with this type of symmetry or 
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greater symmetry. For example a polynomial which is 
fully symmetric is also by definition symmetric in any 
subset of its arguments. A polynomial which is symmet- 
ric in two subsets of its arguments is also symmetric in 
any further sub-division of those subsets. What we re- 
ally require is a generating function for polynomials of 
one particular symmetry only, which we shall denote by 
Ztv,s, and we emphasize that equations ID3l and lD4l fail to 
do this. In appendix IE 5l wc make the statement "greater 
symmetry" more precise; in particular we show that for 
tableaux containing the same number of boxes, a tableau 
A is of "greater symmetry" than a tableau /i if and only 
if S\ > Sp, where S is the corresponding spin eigenvalue 
of the tableau. 

What we call Zjy^s is the generating function for poly- 
nomials which are eigenfunctions of S 2 (note that this 
is independent of S z eigenvalue). Below, Zm,s{<1) rep- 
resents either Z£™ lon (q) defined in equation (|D4[) or 
^jv°s 0n (?) defined in equation (|D3|) . Our derivation of 
the generating function Z^,s is as follows. First, it is 
easy to see that 

which is simply the statement that for the state where 
S = N/2, the wavefunction is fully symmetric (bosons) 
or fully antisymmetric (fcrmions). For the next highest 
state we have: 

Z n (a) = Z n (q) + Z n (q) , 

that is, the polynomials space whose dimension is gen- 
erated by either equation ID3I or equation ID4i is com- 
posed of both fully symmetric polynomials and polyno- 
mials which are symmetric in N — 1 indices. We can 
re-arrange this equation to give an expression for the gen- 
crating function that we want to get at: 

Z n (q) = Z n (q) — Z n (q) ■ 
In general we have 



Z N,S (q) = ^N^-i te)' 

i=0 ' 2 

and so 

Zn,s (q) = Z N , S (q) - Z N , s +i (q) ■ (D5) 

These generating functions are used to calculate the re- 
sults displayed in table HI A more rigorous proof of this 
result is given in appendix IE 51 

Although the generating functions provide us with a 
means to calculate the dimension of the space of polyno- 
mials at a given order, they shed no light on the explicit 
form of the polynomials. An insight into the forms of the 



polynomials can be obtained by first re-writing the gen- 
crating function using the following result (see appendix 
EH): 

N fc max 
ZN,s{q) = q J \{ Y —Y. h ' k 1^ (D6) 

n=2 1 k=l 

where b'k are positive integer coefficients and 

fcmax = (y - S) (y + S) . 

We interpret the expression for the generating func- 
tion in (|D6[) as describing a polynomial basis comprising 
"primitive" polynomials at orders J up to fc max . , each of 
which has b' k linearly independent coefficients. In order 
to indicate the polynomial orders at which the primitive 
polynomials occur, we include the following Y functions 
in table EH] 

v _ z N ,s (q) _ , , k 

Yn,s - — 2^b k q . 

FT 1 k=i 

II l- q l 
1=2 

Some examples of the b'k can thus be read off from table 

rvm 

The conclusion of this section is the following key state- 
ment: the most general polynomial of a given symmetry 
type is a linear combination of the primitive polynomi- 
als multiplied by any valid translationally invariant fully 
symmetric polynomial. 

b. Calculation of Primitive Polynomials 

We shall conclude our argument by describing the pro- 
cedure by which we have calculated the primitive poly- 
nomials presented in tables |TT] and IIIII We start with the 
most general fixed degree translationally invariant poly- 
nomial as the spatial wavefunction: 

iV 

^spatial = ^ a ll ... lN z\ 1 ...z l £ with y^i m =L. 

ii,...i n m— 1 

We then apply the appropriately constructed symmetric 
group algebra projection operator, e^. This will leave a 
sum of terms which arc of the correct permutation sym- 
metry, and hence they can be written in terms of a basis 
of elementary symmetric polynomials: 

$(N,S,1)= bi 1 ...i N e 1 2 1 N (z 1 ...z N )...e l ] $ N (z 1 ...z N ) 

N 

with m hn = L. 

m— 1 

This is equivalent to a basis involving primitive polyno- 
mials and fully symmetric polynomials only. We start 
with linearly independent coefficients a^...^, but once 
we apply the projection operator we have a smaller set 
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TABLE VII: These functions are calculated by taking 
the ratios of the proper generating functions for 
polynomials of symmetry type A to the generating 
function of a fully symmetric polynomial for the same 
number of particles. The polynomial order of each term 
indicates the presence of a primitive polynomial at that 
order. The value of the coefficient indicates the number 
of linearly independent parameters constituting the 
primitive polynomial at that order. 
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of coefficients bi 1 ...i n - These coefficients are not linearly 
independent in general. We can express the b coefficients 
in terms of the a coefficients, and in principle it is possible 
to solve these systems of equations to determine all possi- 
ble linear dependences between the b coefficients. Indeed, 
using this procedure, we have independently verified the 
number of linearly independent coefficients predicted by 
the generating function method. 

If we perform this construction at a polynomial order 
exactly corresponding to a primitive polynomial then we 
can be sure that one of the basis polynomials (i.e., the 
terms multiplying one of the bi t ...i n coefficients) can be 
considered as the primitive polynomial. We can therefore 
extract the primitive polynomial by setting an appropri- 
ate number of the b^,,,^ coefficients to zero, such that 
the remaining expression contains the correct number of 
linearly independent coefficients required for the primi- 
tive polynomial. For N < 5, almost all of the primitive 
polynomials contain only one arbitrary coefficient and so 
this procedure can be reduced to setting all but one of 
the 6j 1 ...i„ coefficients to zero. The polynomials given in 
table HIl arc all derived using this technique. We have elu- 
cidated the details of the technique here in order to point 
out that the primitive polynomials listed in table [XT] arc 
defined in a somewhat arbitrary way. Despite this, once 
we reconstruct the spatial wavefunctions from primitive 
polynomials multiplied by symmetric polynomials with 
a fixed combined degree, then these functions will fully 
span the basis. 



We began this construction procedure with the most 
general possible Spatial- In practice we are free to choose 
a less general starting spatial wavefunction if we know 
what properties to expect in the final spatial wavefunc- 
tion after projection (in other words we can effectively 
pre-empt the linear dependences that will be introduced 
between the coefficients and this saves a great deal of 
computational time). The first part of the spatial wave- 
function will be associated with the first standard tableau 
of the shape A = [Ni, N%, ...Nk], The corresponding poly- 
nomial will be translationally invariant and symmetric 
or antisymmetric in indices 1 to N\, in indices JVi + 1 
to ./V2, and so on. Using this information we can choose 
the starting spatial wavefunction to be the most general 
translationally invariant analytic polynomial with these 
symmetries at a given order, $ re duced- These symmetries 
are in fact uniquely associated with the r = 1 standard 
tableau and as a result it follows that 

reduced- 

(D7) 

We shall now discuss a simple example illustrating 
the whole construction procedure in the case of a two- 
component system. For simplicity there are no undeter- 
mined parameters in this example. Let us consider the 
three-particle L — 1, S — 1/2 , S z — 1/2 fermion wave- 
function. We begin by constructing the spatial wavefunc- 
tion associated with the first spin eigenfunction. To do 
this we first construct the r = 1 standard tableau asso- 
ciated with N ~ 3 and S = 1/2. We know that the spin 
eigenfunction X (3, 1/2, 1/2, 1) is constructed from: 
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3 





The fermion spatial wavefunction associated with the 
spin eigenfunction X(N = 3, S = 1/2, S z = 1/2,1) 
is then constructed from the contragradient standard 
tableaux, which is 
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3 


2 





For bosons we instead construct the spatial wavefunction 
from the uppermost tableau shown here. 

We have already calculated X (3, 1/2, 1/2, 1) in ap- 
pendix IC 11 and it is clearly symmetric under exchange 
of the first two spins, but has no other symmetry. If we 
require that the overall wavefunction is antisymmetric 
(for fermions) then we know immediately that the spa- 
tial wavefunction $(iV = 3, S = 1/2, 1) will at the very 
least be antisymmetric under exchange of the first two 
co-ordinates. Using this information we can choose the 
primitive spatial wavefunction to be 

^reduced = {Z\ - Z 2 ) 

We then apply a projection operator to generate 
$(3, 1/2, 1) (see appendix lEl for the explicit form of this 
operator): 

$(3, 1/2, 1) = eff$ = e [ lf (zt - ~z 2 ) = (zi - h) 
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As we have explained, this term is associated with the 
first primitive spin wavcfunction and we know that all of 
the physics is contained in the construction 

im> (h - h) . 

Given this information we can effectively determine the 
remainder of the wavcfunction by antisymmctrizing in all 
the particle indices (see Eq. |2"2"|) . 

Although simply antisymmctrizing does indeed give 
the correct result, it may be necessary to write the result 
explicitly in terms of spin eigenf unctions. In this exam- 
ple we have / A = 2. We associate $(3, 1/2, 1) with the 
first spin eigenfunction: 

X (3, 1/2, 1/2, 1) = (2 im> - mt) - \ltf)) ■ 

The second spin eigenfunction is: 

A(3,l/2,l/2,2) = ^(mt}-|in». 

The associated spatial wavefunction is given by applying 
the operator (written in terms of permutation operators) 

4 2 i 1] = ^ {(312) - (132) - (231) + (321)} 

to the first part of the spatial wavcfunction. This results 
in the polynomial 

*(3,l/2 J 2) = -(zi + f2-2i 3 )- 

The full wavefunction must be a the sum of these two 
terms: for L = 1 this is 

ii, 1/2, 1/2, i) = (2 mi) - \m) - im» (zi - z 2 ) 
-(\m)-\m))(zi + z2-2z 3 ). 

Comparing with Eq. [22] we see that once the terms cor- 
responding to each of the primitive spin wavefunctions 
are collected up then the results do indeed agree. 

To obtain the equivalent boson wavefunction in this 
case we can simply interchange the associations between 
the spin and spatial parts. We can, of course, apply the 
same technique described in this example to construct the 
wavcfunction at any order, starting from the first part of 
the spatial wavefunction at that order. 

c. Generlization to the Multi- Component Case 

In the multi-component case, we have already shown 
that the full spatial wavcfunction can be reconstructed 
from just the first part, <J>(7V, A, 1). In this section we 
aim to describe the form of this first part of the spatial 
wavcfunction, a translationally invariant analytic poly- 
nomial of a particular symmetry, in terms of the elemen- 
tary symmetric polynomials. We would like to derive 
an expression for the dimensions of the spaces of such 



polynomials, as we did for the two-component case, by 
constructing a set of generating functions. 

We start by writing down a generating function for 
the space of polynomials with an arbitrary number of 
symmetry subsets, and we do this by simply multiply- 
ing together the generating functions corresponding to 
each of the different subsets (previously we only needed 
to consider two subsets). To take into account the trans- 
lational invariance of the polynomial, we use the relative 
co-ordinates also defined in section [TTJ This introduces 
the constraint given in Eq. ID1I To enforce this con- 
straint, we remove one of the sets of e\ terms by starting 
one of the generating functions from a lower bound of 2 
rather than 1 (similar to the argument given in for the 
two-component case). 

N! N 2 N„ 

z boso„ = TT _J_ TT _J_... TT -J— 
A 11 \ (.mi 11 \ — n m 2 11 1 — 0iri2 

mi=2 y m 2 = l m„=l y 

(D8) 

or for conjugate tableaux we would have to take out an 
appropriate Jastrow factor of order: 

n 

J = E 2^ {Nl ~ 1} - 

i=l 

The generating function in this case is: 

Ni N 2 N n 

z^~= q j n — - — n — — ... n — — • 

mi— 2 m2 — l m n — l 

(D9) 

In fact the generating functions written down here de- 
scribe a space of polynomials which has a certain sym- 
metry type, or a greater symmetry. For example, the 
generating function Zx=m\-\ gives the dimension of the 
vector space of polynomials of symmetry type A = [2, 1] 
added to the dimension of the vector space of polyno- 
mials of symmetry type A = [3]. More precisely we say 
that the vector space generated by Z\ oson contains two 
sub vector spaces, each of which corresponds to a partic- 
ular irreducible representation of the symmetric group. 
It follows that for our purposes, determining the vector 
spaces of polynomials associated with irreducible repre- 
sentations of the symmetric group, we are looking for 
the dimensions of these sub vector spaces. If we wish to 
perform this decomposition of vector spaces for any ar- 
bitrary shape of tableau then an important question to 
address is which sub vector spaces are included, in gen- 
eral, and how many times is each vector space included 
(as there is no reason in principle why a sub vector space 
cannot be included multiple times in the decomposition). 
The answer to this question can be found in the under- 
lying mathematical structure of the construction of vec- 
tor spaces corresponding to representations of the sym- 
metric group (objects known as modules in the mathe- 
matical nomenclature); we shall simply state the result 
here, leaving the details of the underlying mathematics 
to appendix IE 51 In order to state the result we must 
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first introduce a notation which makes precise the state- 
ment "greater symmetry"—: when comparing two Young 
tableau shapes // = [/zi, \i2, ■■■/J-n] and A = [Ai, A2, 
we say that /i dominates ( has "greater symmetry" than) 
A, written /i> A if 

A*i + + ■■■ + fJ-i > Ai + A 2 + ... + A, V i > 1, 

with 

/ i i>/ = Ai>„ = 0. 

For example [3, 3] > [2,2, 1, 1]. Note that there are cases 
when tableau are incomparable, such as [3, 3] and [4, 1, 1]. 
One way to avoid such incomparable tableau shapes by 
restricting possible shapes to tableau having no more 
than two rows, which is precisely the two-component 
case. For proof of this result see appendix IE 51 Using 
this notation, the general relation between the Zf oson 
generating functions and the generating functions corre- 
sponding to irreducible representations of the symmetric 
group, Z\, is 

z boson = J2 KpxZp. (D10) 

In this equation the coefficients are called Kostka 
numbers^,. It is clear that equation ID10I is a matrix 
equation. Since the matrix of Kostka numbers is non- 
singular, the relation can be inverted, which gives us our 



key result: 

^ = E(^U^ oson - (on) 

A derivation of this result is given in appendix IE 51 In 
table I Villi we have listed the inverse matrices of Kostka 
numbers for tableau of size 3,4 and 5. Using these num- 
bers, we can construct, for example: 

ry ^boson /7boson ^boson , ^boson 

^[3,1,1] — Z [3,l,l] ~ Z [3,2] - z [4,l] "+" ^[5] 

Notice that if we restrict the allowed tableau shapes to 
have no more than two rows then we recover the special 
case of the result given in equation ID5I An identical 
argument applies for the description of conjugate tableau 

^fcrmion 

We have used these results to calculate the lists of the 
dimensions of vector spaces associated with irreducible 
representations of the symmetric group listed in table U 
for the symmetry types not already included in the SU(2) 
case. 



TABLE VIII: A list of inverse Kostka matrices for different Young tablau sizes. Lists of these matrices up to N = 8 
are given in Ref. HH. 



(a) N=3. 








(b) N= 


=4. 








A = 


[3] [2,1] [1,1,1] 




{K)-l A = 


[4] [3,1] [2,2] [2,1,1] [1,1,1,1] 


M = 


[3] 


1 -1 


L 






1 


-1 





1 


-1 


[2,1] 


1 


2 




[3,1] 





1 


-1 


-1 


2 


[1,1,1] 





L 




[2,2] 








1 


-1 


1 










[2,1,1] 











1 


-3 










[1,1,1,1] 














1 










(c) N=5. 














{K)-{ A = 


[5] [4,1] [3,2] [3,1,1] [2,2,1] [2,1,1,1] [1,1,1,1,1] 






= 


[5] 


1 


-1 


1 




-1 




1 








[4,1] 





1 


-1 -1 1 




1 




-2 








[3,2] 








1 -1 -1 




2 




-2 








[3,1,1] 








1 -1 




-1 




3 








[2,2,1] 








1 




-2 




3 








[2,1,1,1] 













1 




-4 








[1,1,1,1,1] 


















1 
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Appendix E: Further Mathematical Details 

1. Young Operators 

The concept of Young operators may be familiar to 
many readers. The Young operators are idempotent pro- 
jection operators which generate minimal left ideals in 
the regular representation of the symmetric group; in 
other words the Young operators form a resolution of 
the identity which generates the irreducible representa- 
tions of the symmetric group algebra^. By definition 
this means that they form a basis spanning the group 
algebra. 

The Young operators are constructed as follows: with 
each standard tableau, T r , we associate an operator 

Here P r is product over all rows in the tableau of the 
sum of all permutations that permute the numbers in 
the same row. Q r is product over all columns in the 
tableau of the sum of all permutations that permute the 
numbers in the same column multiplied by the sign of 
those permutations. Essentially, P r symmetrizes in the 
row indices and Q r antisymmetrizes in the column in- 
dices. The Young operator is the product of these two 
operators. 

The set of Young operators can be used to generate all 
of the results derived in appendix [B] however the result- 
ing projected functions do not form a fully orthogonal 
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It is more convenient therefore to use a set of 



operators which automatically project out an orthogonal 
basis of generalized spin wavef unctions. 

The advantage of the matric basis, defined in appendix 
lAl is that it enables us to use the representation matri- 
ces of any irreducible representation to construct a corre- 
sponding basis of the group algebra, and for convenience 
we choose the orthogonal representation to ensure or- 
thogonal basis functions. 



2. Matric Basis Operators 

In table IIXI below we list the operators e*i that are 
applied to the first part of the spatial wavefunction in 
order to generate the spatial part of the wavefunction as- 
sociated with the r th spin eigenfunction. The operators 
are written in terms of permutation operators: we use 
the notation explained in Ref. HH, for example. In this 
notation the term (p\,P2, ■ ■■,Pn) denotes the permuta- 
tion operation that relabels index 1 with index pi , index 
2 with p2 and so on. 

The operators here are calculated using a 
Mathematica™-based computer program^. This 
program can be used to calculate expressions for 
operators up to N = 5 and for any valid A. 



In appendix IC 41 we discussed the construction of spin 
eigenfunctions describing systems containing two types of 
electrons corresponding to two valleys labeled by 1 and 
2 or A and B. 

The general mathematical procedure for the construc- 
tion of spin eigenfunctions corresponding to multiple sets 
of spinful particles is explained in detail in Ref. l68l . The 
procedure provides a group-theoretical explanation of the 
results derived by the brute-force approach presented in 
appendix IC 41 What we are doing is starting with a 
decomposition of wavefunctions into representations of 
SU(n) and then further decomposing into representations 
of SU(2) x SU(2) (for two sets). In general we could 
in principle decompose into representations of SU(a) x 
SU(&) x ... . That would describe a system of a compo- 
nent particles in set 1, a system of b component particles 
in set 2 and so on. 

Let us illustrate the method with the following ex- 
ample: decomposing a four-component system into two 
two-component systems (which applies to the case of 
graphenc when the system is spin rotationally invariant). 
Given a number of electrons, N, in the system we parti- 
tion this number into two parts (representing the number 
of electrons in each set). For example, with N = 3 we 
can have the partitions [3,0], [2,1], [1,2] and [0,3]. Next 
we construct SU(2) spin eigenfunctions corresponding to 
each set of partitions. In this example we would con- 
struct spin eigenfunctions corresponding to the following 
SU(2) Young tableaux shapes: 



[3,0] 



5 = 3/2 ; 5 = 



[2,1] 



5=1/2 ; 5 = 



□ 



5 = 1 ; 5 = 1/2 



; □ 



5 = ; 5 = 1/2 

and vice versa for the [0,3] and [1,2] partitions. The cor- 
responding SU(2) spin angular momentum eigenvalues 
are written below the Young tableaux. The next step 
is to construct SU(4) Young tableaux by forming tensor 
products between all combinations of the SU(2) Young 
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tableaux for each partition. This is done accordin g to 
the Littlewood Richardson rule, see for example Ref. |48|. 
Simultaneously we add the angular momenta of the two 
sets for each partition. For the [3,0] and [0,3] partitions 
this is a trivial procedure since there is only one set con- 
taining electrons. For the [2,1] we construct: 



4. Generating Functions for Symmetric 
Polynomials 

Properties of Generating Functions 

The generating function for a translationally invariant 
symmetric polynomial in N variables is given by 



-®n= 





— 1© 







N 1 

z N {q)= n — 



and we add the angular momenta according to the stan- 
dard procedure: + 1/2 = 1/2. Also we have: 



Consider a product of two generating functions for sym- 
metric polynomials 

Z Nl (q) Z N2 (q) . 



and 1 + 1/2 = 3/2, 1/2. The final step is to associate all 
possible shapes of SU(4) tableaux deduced from combin- 
ing each pair of cluster tableaux with all possible results 
of the corresponding addition of angular momenta. For 
example, the SU(4) tableaux shape A = [2, 1] can occur 
with angular momenta S = 3/2 in 2 places (once in clus- 
ter partition [2,1] and once in partition [1,2] ). The same 
shape can occur with angular momentum S = 1/2 in 6 
places (once in [3,0] and [0,3] and twice in [2,1] and [1,2]). 
This result matches our previous calculation: the shape 
A = [3] occurs with S = 3/2 in 4 places and with 5* = 1/2 
in 2 places. Finally, the shape A = [1,1,1] occurs with 
only with S = 1/2 in 2 places. We have checked that all 
of these results are consistent with the alternative brute- 
force method described in appendix I C 41 



In general one might want to consider a system in 
which there are more than two sets of particles. We 
might also wish to consider sets of spin-j particles. In 
principle the procedure for constructing the basis of psc- 
duopotentials is a simple generalisation of the procedure 
described here for the case of graphene: we write down 
all valid symmetry types for each set of cluster partitions 
and all the corresponding angular momentum eigenval- 
ues. For example, if we have 3 sets of electrons, then we 
would write down the set of all SU(2) tableaux for each 
partition of N into the three sets, and then construct all 
possible combinations of tensor products between these 
tableaux. In the case where j ^ 1/2 some additional 
complications arise due to the multiple angular momen- 
tum eigenvalues associated with each tableaux shape (the 
result of what is called fractional parentage). The reso- 
lution of these issues is described in Ref. K38l 



The dimension of the space of polynomials at order L is 
given by: 



d(L,N u N 2 ) 



1 

ZJ \dq 



Z Nl (q) Z N2 (q) 



9=0 



but we can expand this expression out using the Leibniz 
rule: 



z~2 L ° k 



k=0 



L-k 



-'JVi- 



So the dimension of the space of the combined symmetry 
is given by the sum of the dimensions of the spaces of 
polynomials of the same symmetry but with all possible 
subdivisions of the angular momentum, L. This is pre- 
cisely the result we require for the generating function of 
a polynomial with two combined symmetries. This re- 
sult clearly generalizes to a product of arbitrarily many 
generating functions. 

Now consider a modified generating function given by 
multiplying by a factor of q J . The dimension of the 
space of polynomials arising from this modified gener- 
ating function is effectively a special case of the above 
result: 



d(L,N) 



fc=0 



(*) 



L-k / \ k 
q j( d 



tfC L -jJ\[± 



(L-J)! \dq 



dq 
L-J 



L-J 



m=2 
N 

n ^ 

rn=2 



dq 
N 

n Th 



9=0 
N 

n 

m=2 



9=0 



9=0 



d sym (L - J, N) 



9=0 
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Dividing Generating Functions 
In appendix [Dl we made use of the following results: 



Z N ,s(q) 

N 

g J Tl 

n = 2 



+ S J$—S 

n irjr n i=V 



I 

/ = 2 



n 

n = 2 

V (i-q w+1 - J ) 

3=1 



A< 



at 



S 



The final expression is a so-called q-binomial 
coefficient^. We can thus simplify the generating func- 
tion to 

N 1 

n=2 q k 

where bk are positive integer coefficients. Using the defin- 
tion of q-binomials we can show 



Using the result from equation (|D5|) we have 



Zn, s (q) = j 

N q 

n 

n=2 

It follows that 



JV 



N 

-S 



,1 + J+2S 



N 
-5-1 



71=2 



fe=l 



where b'k are positive integer coefficients, which can be 
derived from the previous formula above in terms of the 
q-binomial coefficients, and 



S) (f + S) 



5. Vector Spaces Associated with Representations 
of the Symmetric Group 

In this final appendix we shall present an alternative 
derivation of the results presented in appendix [D] Our 
aim is to construct a basis of polynomials which are in 
a one-to-one correspondence with irreducible representa- 
tions of the symmetric group, and hence are labeled by 
Young tableau shapes A = [JVi,iV 2j ...N n ]. The starting 
point of our derivation is the set of symmetric polynomi- 
als, which have subsets of symmetries. Such polynomi- 
als can also be enumerated in terms of integer partitions 
and hence Young tableau shapes. The question is how 
are these polynomial constructions related. In this ap- 
pendix we shall explain how this question can be posed 
in a more precise mathematical sense, in terms of the 



representation theory of the symmetric group. Having 
demonstrated an equivalence to a problem in mathemat- 
ics, we shall describe the solution and then finally explain 
how this solution can be applied in the context of our 
problem. 

To begin with, recall the simple argument given in ap- 
pendix [Dl which applied to the case of Young tableaux 
restricted to having no more than two rows. In that 
case we stated that the basis of polynomials of a given 
symmetry type actually contains polynomials of that 
particular symmetry type or of a greater symmetry 
type. We should make this statement more precise. 
Mathematically speaking the term "greater symmetry" 
is equivalent to whether one tableau shape dominates 
another—. When comparing two Young tableau shapes 
fi = \fJLi,fj,2,...iJ,n\ and A = [Ai, A 2 , ...A;], we say that 
fj, dominates ( has "greater symmetry" than) A, written 
fi>\ if 



Mi + M2 + ••• + Mi > Ai + A 2 + ... + Aj V i > 1, 



with 



Mi>! — Ai>„ — 0. 



In the special case where the Young tableau A containins 
N boxes in two rows, with N/2 + S\ boxes in the first 
row, this definition reduces to just: 



Mi > Ai => f 
Mi + M2 > Ai + A 2 



N 



— 2 



N> N. 



The second statement is clearly always satisfied, and the 
first statement is true for all spin eigenvalues S\ < S^, 
where S\ is the spin eigenvalues associated with sym- 
metry type A. This is precisely the way in which we 
defined "greater symmetry" in appendix [D] An impor- 
tant consequence of this result is that there are no pos- 
sible cases of ambiguous dominance (that is, where two 
tableau shapes do satisfy some but not all of the condi- 
tions for dominance, for example [3,3] and [4,1,1]). We 
can represent the set of dominances diagrammatically in 
a Hassc diagrarr>2U. For the case of Young tableau con- 
taining TV boxes in only two rows the Hasse diagram is: 

[TV] 
I 

[AT -1,1] 
I 

[N - 2, 2] 

I 

[N/2 + S^,N/2-S^] 

I 

I 

[1,1,1,...] 

The argument given in appendix [D] to derive equation 
(|D5jl follows directly form the ordered structure of this 
set of dominaces. 
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More generally, if we want to consider arbitrary shapes 
of Young tableau, the set of dominances becomes more 
complicated. For example the Hasse diagram for all 
N = 5 Young tableau is: 

[5] 

I 
[4,1] 

I 
[3,2] 

I 

[3,1,1] 

S \ 
[2,2,1] [2,1,1,1] 

\ S 
[1,1,1,1,1] 

In general we must be more careful take into account 
the exact structure of the inclusions each of these sym- 
metry subsets, and in particular we cannot rule out that 
in general one subset may be included multiple times 
within a subset of greater symmetry. At this juncture we 
shall introduce the equivalent mathematical formulation 
of this problem. 

In the mathematical theory of representations of the 
symmetric group the symbol M x denotes the vector space 
spanned by objects which correspond to the set of all 
Young tableau of shape A. For example M^' 1 ^ has the 
basis 



2 


3 | 4 | 


1 


3 | 4 | 


1 


2 | 4 | 


1 


2 | 3 | 










_3_ 




4_ 





We say that M forms a representation of the symmet- 
ric group called a permutation module. In the context 
of polynomials, the permutation module corresponds ex- 
actly to the vector space of polynomials associated with 
a particular symmetry type. This is because the symme- 
try type of a polynomial should correspond to a Young 
tableau shape, as we have already explained, and because 
there should be no restriction on the order of the indices 
in a symmetric polynomial, so we must take into account 
all possible ways to place the numbers in a tableau of 
a given shape. The dimension of this polynomial vector 



space is thus given by the generating functions which we 
have denoted by Z\. 

Our aim is to construct the space of polynomials corre- 
sponding to irreducible representations of the symmetric 
group. In the mathematical theory of representations of 
the symmetric group the symbol 3 denotes the vector 
spaces which correspond to irreducible representations of 
the symmetric group (the notation is not to be confused 
with our S\, which is the spin eigenvalue associated with 
symmetry type A). These vector spaces are spanned by 
objects which correspond to the standard Young tableau 
of shape A. For example S^ 3 ' 1 ^ has the basis 



1 


3 | 4 | 


1 


2 | 4 | 


1 


2 | 3 | 


■2_ 




J3_ 









S x are called Specht modules^. In the context of poly- 
nomials, these Specht modules correspond exactly to the 
space of polynomials that becomes associated with irre- 
ducible representations of the symmetric group, in other 
words, precisely the space of polynomials that we require 
for our investigation. The dimensions of this space of 
polynomials are given by the generating functions which 
we have denoted by Z\. 

A well studied mathematical problem is to decompose 
the permutation module M x in terms of irreducible rep- 
resentations of the symmetric group, or equivalently, in 
terms of Specht modules. It turns out that it is possible 
to write M A as a direct sum of Specht modules, and the 
result isZS: 

M x = © K^xS*. 

The coefficients appearing here are called Kostka 
numbers^. 

Finally, since this relation involves a direct sum of ir- 
reducible representations, it follows that the dimensions 
of these irreducible representations simply add together, 
and we find: 

It also follows that an equivalent relation must also be 
satisfied by the corresponding generating functions, and 
hence we find equation (jDlOp . 
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TABLE IX: Some examples of the matric unit operators. These operators form an orthogonal basis of the symmetric 
group algebra in correspondence with the orthogonal representation of the symmetric group. 

(a) Matric unit operators associated with the orthogonal representation of the symmetric group. 



en 11 = H 2 ( 123 ) + 2 ( 213 ) ~ ( 312 ) - ( 132 ) - ( 231 ) ~ ( 321 )> 

= {( 312 ) + ( 132 ) - ( 231 ) - ( 321 )> 

{ 3 (1234) + 3 (2134) + 3 (3124) + 3 (1324) + 3 (2314) + 3 (3214) - (4213) - 
(2413) - (1423) - (4123) - (2143) - (1243) - (1342) - (3142) - (4132) 
- (1432) - (3412) - (4312) - (4321) - (3421) - (2431) - (4231) - (3241) - (2341) _ 

-[3,1] _ V2 / 2 ( 4213 ) + 2 ( 2413 ) + 2 ( 1423 ) + 2 ( 4123 ) + 2 ( 2143 ) + 2 ( 1243 ) ~ ( 1342 ) - ( 3142 ) - ( 4132 ) 
621 3 1 - (1432) - (3412) - (4312) - (4321) - (3421) - (2431) - (4231) - (3241) - (2341) 




(1342) + (3142) + (4132) + (1432) + (3412) + (4312) 1 
- (4321) - (3421) - (2431) - (4231) - (3241) - (2341) J 



{ 2 (1234) + 2 (2134) - (3124) - (1324) - (2314) - (3214) - (4213) - (2413) 
- (1423) - (4123) + 2 (2143) + 2 (1243) - (1342) - (3142) - (4132) - (1432) 
+2 (3412) + 2 (4312) + 2 (4321) + 2 (3421) - (2431) - (4231) - (3241) - (2341) 

,[2,2] _ 73 / ( 3124 ) + (1324) - (2314) - (3214) + (4213) + (2413) - (1423) - (4123) 1 
621 ~~ ~ | + (1342) + (3142) - (4132) - (1432) + (2431) + (4231) - (3241) - (2341) J 



(b) Matric unit operators associated with the contragradient orthogonal representation of the sym- 
metric group. 



4 2 ' 11 = i {2 (123) - 2 (213) - (312) + (132) - (231) + (321)} 
4i 1] = # {( 312 ) - ( 132 ) - ( 231 ) + ( 321 )> 



*[3,1] 

t-1 1 



3 (1234) - 3 (2134) + 3 (3124) - 3 (1324) + 3 (2314) - 3 (3214) - (4213) + 
v (2413) - (1423) + (4123) - (2143) + (1243) - (1342) + (3142) - (4132) 

+ (1432) - (3412) + (4312) - (4321) + (3421) - (2431) + (4231) - (3241) + (2341) 



^5 I 2 (4213) - 2 (2413) + 2 (1423) - 2 (4123) + 2 (2143) - 2 (1243) - (1342) + (3142) - (4132) 
~ ' + (1432) - (3412) + (4312) - (4321) + (3421) - (2431) + (4231) - (3241) + (2341) 



3 [3,1] 



(1342) - (3142) + (4132) - (1432) + (3412) - (4312) 
- (4321) + (3421) - (2431) + (4231) - (3241) + (2341) 



{ 2 (1234) - 2 (2134) - (3124) + (1324) - (2314) + (3214) - (4213) + (2413) 
- (1423) + (4123) + 2 (2143) - 2 (1243) - (1342) + (3142) - (4132) + (1432) 
+2 (3412) - 2 (4312) + 2 (4321) - 2 (3421) - (2431) + (4231) - (3241) + (2341) 

~[275] _ v| f (3124) - (1324) - (2314) + (3214) + (4213) - (2413) - (1423) + (4123) 1 
621 2 I + (1342) - (3142) - (4132) + (1432) + (2431) - (4231) - (3241) + (2341) J 
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